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ABSTRACT 

We have used the Very Large Array to image the Hi 21-cm hne emission in 
the edge-on Sd galaxy IC 2233 and the blue compact dwarf NGC 2537. We also 
present new optical B, R, and Ha imaging of IC 2233 obtained with the WIYN 
telescope. Despite evidence of localized massive star formation in the form of 
prominent H II regions and shells, supergiant stars, and a blue integrated color, 
IC 2233 has a low surface brightness disk with a very low global star formation 
rate ( <O.O5M0 yr~^), and no significant 21-cm radio continuum emission. The 
H I and ionized gas disks of IC 2233 are clumpy and vertically distended, with 
scale heights comparable to that of the young stellar disk. Both the stellar and H I 
disks of IC 2233 appear flared, and we also find a vertically extended, rotationally 
anomalous component of Hi extending to ~2.4 diQ kpc from the midplane. The 
H I disk exhibits a mild lopsidedness as well as a global corrugation pattern with 
a period of ~ 7 dio kpc and an amplitude of ~ 150 dio pc. To our knowledge, this 
is the first time corrugations of the gas disk have been reported in an external 
galaxy; these undulations may be linked to bending instabilities or to underlying 
spiral structure and suggest that the disk is largely self-gravitating. 

Lying at a projected distance of 16'7 from IC 2233, NGC 2537 has an Hi disk 
with a bright, tilted inner ring and a flocculent, dynamically cold outer region that 
extends to ~3.5 times the extent of the stellar light (-D25)- Although NGC 2537 
is rotationally-dominated, we measure H I velocity dispersions as high as o"v,hi ~ 
25 km s~^ near its center, indicative of significant turbulent motions. The inner 
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rotation curve rises steeply, implying a strong central mass concentration. Our 
data indicate that IC 2233 and NGC 2537 do not constitute a bound pair and 
most likely lie at different distances. We also find no compelling evidence of a 
recent minor merger in either IC 2233 or NGC 2537, suggesting that both are 
examples of small disk galaxies evolving in relative isolation. 

Subject headings: galaxies: spiral — galaxies: fundamental parameters — galaxies: 
kinematics and dynamics — galaxies: individual (IC 2233, NGC 2537) — galaxies: 
ISM— ISM: Hi 



1. Introduction 

This is the second in our series of papers presenting H I imaging observations of edge- 
on, pure-disk galaxies obtained with the Very Large Array (VLA0 We are focusing on 
'superthin' Sd galaxies that are characterized by their highly flattened, dynamically cold 
stellar disks and the absence of a visible spheroid component (Goad & Roberts 1981). The 
structural simplicity and edge-on orientation of superthins makes these galaxies excellent 
laboratories for exploring how internal versus external processes regulate many of the key 
properties of late-type disk galaxies, including their structures and star formation histories. 
For example, the abundance of thin, bulgeless disks in the local universe (Karachentsev et 
al. 1999; Kautsch et al. 2006) poses challenges for hierarchical galaxy-formation models in 
which galaxies are built-up through violent mergers (Abadi et al. 2003; D'Onghia & Burk- 
ert 2004; Eliche-Moral et al. 2006; Kormendy & Fisher 2005 and references therein). The 
superthin disks that we are studying also inhabit a particularly interesting galactic mass 
range iV^ot ~100 km s~^), where important changes in the structure and composition of the 
interstellar medium have been noted (Dalcanton et al. 2004; Matthews et al. 2005). Com- 
plementary studies of both stellar and interstellar components are necessary to understand 
the physics behind these changes and, in turn, how they affect the regulation of star forma- 
tion. Finally, it is well established that the thinnest galaxies tend to be chemically unevolved 
systems of low optical surface brightness (Goad & Roberts 1981; Bergvall & Ronnback 1995; 
Matthews et al. 1999; Gerritsen & de Blok 1999). Thus, the vertical structure of superthin 
galaxies provides information relevant to understanding the overall disk structure and mass 
distribution of low surface brightness (LSB) galaxies in general. 
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In a previous paper (Uson & Matthews 2003; hereafter Paper I), we presented VLA 
observations of the isolated, LSB superthin spiral UGC 7321. Here we turn to IC 2233 (= 
UGC 4278). In addition to new VLA Hi observations we present optical B, R and 
Ha imaging of IC 2233 obtained with the WIYN telescope (§3]). IC 2233 has a comparable 
physical size and luminosity to UGC 7321, and both are weak infrared sources with very low 
current star formation rates. However, IC 2233 is slightly less massive and has a slightly 
higher optical surface brightness. And while both systems are H l-rich and exhibit bulgeless, 
highly flattened stellar disks with minimal dust absorption, our new optical and H I imaging 
study has uncovered some interesting differences in the kinematics of these two galaxies and 
in the structures of their ionized and neutral gas disks. 

One of the goals of our VLA study of superthin, edge-on spiral galaxies is to explore the 
role of environment in shaping the structure, kinematics, and evolution of these dynamically 
cold and seemingly "fragile" disk systems. While UGC 7321 appears to be an extremely 
isolated galaxy (Paper I), we targeted IC 2233 because of the presence of a second galaxy at 
a projected separation of 16'7 and a radial velocity difference of AV^ ^ 114 km s~^: the blue 
compact dwarf galaxy (BCD) NGC 2537 (= UGC 4274 = Arp 6 = Mrk 86 = "The Bear's 
Paw"). We observed NGC 2537 simultaneously with the VLA, and our new observations 
highlight several interesting features of this galaxy ( §8.2p . However, we find that IC 2233 
and NGC 2537 do not appear to be physically associated and that a mutual tidal interaction 
is unlikely to have shaped the properties of either galaxy. Indeed, like UGC 7321, both 
IC 2233 and NGC 2537 appear to be examples of low-mass disk galaxies evolving in relative 
isolation f^ TOlTOD . 



2. The IC 2233 and NGC 2537 Field 

The existence of a true physical association between IC 2233 and NGC 2537 has been 
a matter of previous debate (see de Vaucouleurs et al. 1976; Schneider & Salpeter 1992; 
Oosterloo 1993; Gil de Paz et al. 2000a, Wilcots & Prescott 2004). Their surrounding field 
also contains a number of fainter galaxies whose relationship to our targets has remained 
unclear. Recently, some new insight has been provided by observations from the Sloan Digital 
Sky Survey (SDSS)@ 
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We have obtained the pubhcally available spectra and five-band photometry of the 
galaxies observed by the SDSS within one degree of the position of IC 2233. Figure [1] 
is a grey-scale representation of a composite u, g, r mosaic frame of the field with the 
SDSS redshifts of the galaxies indicated. The SDSS spectra confirm the redshifts of IC 2233 
(0.0019) and NGC 2537 (0.0015), but also reveal a group of galaxies at 2; ~ 0.04 that includes 
NGC 2537A, whose association with NGC 2537 had been uncertain (see Davis & Seaquist 
1983, Gil de Paz et al. 2000a). It now appears likely that a prior report of a redshift for 
NGC 2537A equal to that of NGC 2537 was due to confusion of the two galaxies in a single- 
dish Hi measurement. Indeed, we have not detected any Hi emission toward NGC 2537A 
with the VLA. As described in §9.H we have also performed a systematic search for other 
gas-rich neighbors to IC 2233 and NGC 2537, but have not identified any candidates. 

The association between IC 2233 and NGC 2537 has remained controversial because 
previous distance estimates for both galaxies span a range of values. Assuming that the 
galaxies are at rest with respect to the cosmic background radiation leads to corrected 
recessional velocities of ~ 710 km s^\lC 2233) and ~ 590 km s'^NGC 2537) and Hubble- 
fiow distances of ~ 10.0 Mpc (IC 2233) and ~ 8.3 Mpc (NGC 2537), with a formal uncertainty 
of 5% (assuming a Hubble constant of 71 km s^^ Mpc~^). 

Estimates of the distance to IC 2233 derived from velocity fiow models range from 10.4- 
10.9 Mpc (Tully 1988; Schwarzkopf & Dettmar 2000; Garcia- Ruiz et al. 2002). From the 
5-band Tully-Fisher relation, Garcia-Ruiz et al. (2002) derived a distance of 8.1 Mpc, while, 
using the /-band Tully-Fisher relation, Gil de Paz et al. (2000a) stated that IC 2233 is located 
at a distance of 13 Mpc (quoting an unpublished, private communication; these estimates 
all assumed a Hubble constant of 75 km s~^ Mpc~^). Based on Hubble Space Telescope 
(HST) data from Seth et al. (2005), Tikhonov et al. (2006) recently derived a distance to 
IC 2233 of 10.4 ± 0.4 Mpc using the "tip of the red-giant branch" method. However, the 
uncertainty quoted by these latter authors appears to be purely statistical, without allowing 
for the uncertainties in the zero-point calibration and metallicity corrections. Re-analyzing 
the same data, A. Seth (private communication) finds a distance modulus of 30.25 using an 
edge-detection algorithm and 30.11 using a Monte-Carlo approach (see Seth et al. 2005 for 
details on these methods). The difference is indicative of the low signal-to- noise of the data. 
Taking into account the uncertainty in the zero-point calibration and metallicity corrections 
therefore raises the uncertainty in the distance to IC 2233 determined with this method to 
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~ 1.5 Mpc. 

In the case of NGC 2537, Tully (1988) estimated a distance of 9 Mpc based on a peculiar 
velocity flow model of the local supercluster. More recently, Sharina et al. (1999) derived a 
distance of 6.9ti;2 Mpc based on the magnitudes of its brightest resolved stars. 

Given the uncertainties in the distances to both NGC 2537 and IC 2233, we cannot 
exclude with absolute certainty the possibility of a true physical association between them. 
However, as we argue in §9.11 even if they do lie at the same distance, it appears unlikely that 
either galaxy has been significantly influenced by a mutual interaction. For the present work 
we adopt d =10 Mpc for IC 2233 and d=7 Mpc for NGC 2537, but throughout the paper, 
we indicate the distance dependence of all physical quantities using the symbols dio and d^ 
to denote the true distances to the galaxies in units of 10 Mpc and 7 Mpc, respectively. 



3. Optical Imaging and Photometry of IC 2233 

3.1. Observations and Data Reduction 

Harris B and R CCD imaging observations of IC 2233 were obtained on 1997 Novem- 
ber 7 using the 3.5-m WIYN telescop^ at Kitt Peak, AZ. Deeper B and R images, as well 
as a narrow-band Hq;-|-[Nii] image were obtained with the same observing setup on 1999 
January 19 under non-photometric conditions. The Hq;+[Nii] image was obtained with the 
WIYN W^015 filter with a bandpass of 73A, centered at 6570A. 

The imaging camera employed a thinned STIS 2048x2048 CCD with 0".2 pixels, giving 
a field of view of ~6'8 per side. The gain was 2.8e~ per ADU and the readout noise was 
~ 8e~ per pixel. Exposure times were 550 seconds in B and 180 seconds in R during the 
1997 observations, and 900 seconds in B, 750 seconds in R, and 1200 seconds in Ha during 
the 1999 observations. Seeing during the 1997 run was 0"7 and during the 1999 run was 
0"5-0"7. 

The images were processed using standard IRAfQ routines. Overscan and bias levels 
were subtracted from individual frames and the data were fiat-fielded using an average of 5 
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dome flats taken with the appropriate filter. This yielded images flat to ~ 2% in B and to 
better than 1% in R. A single bad column in the CCD was corrected via interpolation from 
the adjacent two columns. 

During the 1997 November run we derived calibrations by observing standard fields from 
Landolt (1992) at three different air masses. The data show excess scatter with respect to 
the best-fit secant law, which suggests that some cirrus might have been present during part 
of the night. These uncertainties are reflected in the rather large errors in the photometry 
listed in Table 1. 

Photometry of IC 2233 was performed using a single elliptical aperture. Foreground 
stars, cosmic rays, and background galaxies were removed first via background interpolation 
from surrounding regions. Sky brightness levels were determined by measuring the mean sky 
counts in several boxes at different locations on the image. 

Basic photometric parameters for IC 2233 are presented in Table 1. Total errors in 
the magnitude measurements were computed following Matthews & Gallagher (1997) and 
take into account sky, fiat-field, and Poisson errors as well as the formal uncertainty in 
the photometric solution. They are dominated by the calibration uncertainties described in 
the previous paragraph. We have corrected all quantities for Galactic extinction following 
Schlegel et al. (1998). Our new measurements (excluding corrections for Galactic extinc- 
tion) agree within the uncertainties with those pubhshed by de Vaucouleurs et al. (1991; 
mB=13.48±0.18) and Swaters & Balcells (2002; mij=12.53±0.02). 
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Table 1. Optical and Infrared Properties of IC 2233 



Parameter 


Value 


Ref. 


a (J2000.0) 
<5 (J2000.0) 
Hubble type 
Distance*' 

Ab (mag) 
An (mag) 


08 13 58.9 
+45 44 34.3 
Sd 
8-12 Mpc 
0.223 
0.138 


1 
1 

2 
2 
3 
3 


Measured Quantities 


Position angle 
Inclination 

■mB"''^ 
S - 


172° ± 1° 
88.5° ± 1.5° 
7.0 
5.'17 
13.05 ±0.14 
12.38 ±0.07 
0.67 ±0.15 
21.3 mag arcsec"^ 
32"± 2" 
29" ± 2" 
4'.'9± O'.'l 


2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 


Derived Quantities 


AR,25.i^ 

Mb' 

Lb' 

Lhu 

Lfir 


15.0 dio kpc 
-17.65 

l.SxlO^'d^ij L© 
~ 22.6 mag arcsec^^ 
6.1 ± 1.8 X 10=^''dfo erg s^^ 
7.1 ±0.9 X lO^d^o Lq 


2 
2 

2 
2 
2 

1,2 



Note. — Units of right ascension are hours, minutes, 
and seconds, and units of declination are degrees, ar- 

cminutes, and arcseconds. 

''We adopt a distance of 10 Mpc (see text). Distance- 
dependent quantities are scaled in terms of d\Q , the ac- 
tual distance in units of 10 Mpc. 

''Measured at an ij-band surface brightness of 25.5 
mag arcsec"^. 

°Error estimate includes ±0.08 mag formal error and 
0.11 mag zero point uncertainty. 

''Corrected for foreground extinction only. 

''Error estimate includes ±0.04 mag formal error and 
0.05 mag zero point uncertainty. 

^Measured central surface brightness corrected for 
foreground extinction 

^Disk scale length based on truncated exponential fit. 
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'^Disk scale height based on exponential fit along the 
minor axis. 

'Corrected for internal extinction (0.70 mag) follow- 
ing TuUy et al. 1998. 

''Deprojected central surface brightness based on 
truncated exponential fit including finite disk thickness 
and corrected for internal and foreground extinction. 

'Uncorrected for the contribution of [Nil]. 

References. — (1) NED database; (2) this work; (3) 
Schlegel et al. 1998. 
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3.2. Properties of the IC 2233 Stellar Disk 

Our i?-band image of IC 2233 is presented in Figure [2] (left panel). Here we draw 
attention to a few key features of this image. 

Examining first the central regions of IC 2233, we see that tlie galaxy sliows no well- 
defined equatorial dust lane. Only a few small patches of extinction are evident in projection 
along the inner few kiloparsecs of the galaxy, together with a narrow dust structure inclined 
roughly 45° to the disk plane. This latter feature can be traced both above and below the 
mid-plane (Figure [2]), and its center corresponds with the kinematic center of the galaxy (see 
§6.1.ip . The nature of this minor axis dust lane is presently unclear; it could be dusty debris 
recently accreted onto IC 2233, or it may comprise material being expelled from a starburst 
nucleus or AGN, analogous to what is observed in M82 (Alton et al. 1999). However, Seth 
et al. (2006) reported no evidence for a nucleus in IC 2233. Another possibility is that it 
represents dense gas streaming along a vertical density perturbation such as a spiral arm. 

Because of the lack of a true equatorial dust lane in IC 2233, it is difficult to determine 
precisely the inclination of the galaxy. Using the 25.5 mag arcsec"^ i?-band isophote (after 
correction for foreground extinction), we measure a semi- minor to semi- major axis ratio 
6/a=0.14. This ratio is very close to the mean intrinsic ffattening expected for Sd galaxies 
[(6/a)o ~0.15; Yuan & Zhu 2004], suggesting that IC 2233 is observed no more than a degree 
or two from edge-on. This is consistent with the small scale height of the stellar disk (see 
below) and the symmetry in the light distribution about the midplane. We therefore adopt 
i =88°.5±1°.5. 

We have estimated the radial scale length of the IC 2233 disk by fitting a truncated 
exponential model to the major-axis light profiles in the B and R bands. In the i?-band 
we derive hr^B = 32" ± 2" (~1.6 dio kpc) with a cutoff radiu^ of 165" and in the R- 
band we find = 29" ± 2" (~1.4 diQ kpc) with a cutoff radius of 185". We caution 
however that an exponential does not provide an ideal fit to the major axis light profile of 
IC 2233; it significantly underestimates the amount of light in the intervals r = 0" — 10" and 
r = 90" — 140" in B and r = 65" — 140" in R. Our scale length estimates should therefore be 
considered only as approximations to the more complex underlying stellar light distribution. 

Perpendicular to the disk, the i?-band brightness profile along the minor axis of IC 2233 
is well described by an exponential function with scale height 4"9±0"1 (240 dio pc). This is 
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somewhat larger than the i^'s-band scale height hz^K='i"-77 reported by Seth et al. (2005B 
However, this difference is consistent with other edge-on, late-type galaxies where scale 
heights are found to be systematically larger in i?-band compared with the near-infrared 
(Bizyaev & Mitronova 2002). In agreement with its identification as a superthin galaxy by 
Goad & Roberts (1981), the global scale height of IC 2233 is on the small end for Sd galaxies 
(cf. Bizyaev & Mitronova 2002; Seth et al. 2005), although it is notably thicker than the 
more massive superthin galaxy UGC 7321 {h^ ~ 150 diQ pc; Matthews 2000). In addition, we 
see evidence that the stellar disk of IC 2233 flares at larger galactocentric radii (see below). 

The excellent seeing during our 1999 observing run shows the disk of IC 2233 par- 
tially resolved into a number of bright, point-like sources. These features are most readily 
identifiable in the outer reaches of the disk, where they stand out strongly against the com- 
paratively faint, diffuse underlying disk light. Our Ha image (discussed below) does not 
reveal any compact sources at these locations, implying that the majority of these are not 
compact H ll regions. Using aperture photometry, we have estimated magnitudes and colors 
for 38 point sources throughout the galaxy that appear to be unblended and do not cor- 
respond with bright, underlying Ha emission. After correction for foreground and internal 
extinction, the colors we estimate for these features lie in the range {B — R)q 0.2 to 3.5, and 
absolute magnitudes range from Mr ^ —7 to —8.7. These ranges of color and luminosity are 
consistent with a population of blue and red supergiant stars (Johnson 1966; Dohm-Palmer 
& Skillman 2002). While some of the bluer sources may be young or intermediate-age star 
clusters or associations, 15 of the sources have {B — R)q > 2, making them redder than 
even the oldest globular clusters (e.g.. Bridges et al. 1996). This is another indication that 
we are resolving a population of individual luminous stars. Using HST images, Seth et al. 
(2005) also identified a population of red and blue supergiants in a portion of the inner disk 
of IC 2233. 

Supergiants represent massive stars that have depleted their core hydrogen, and those 
luminous enough to be identified in our ground-based images are likely to have ages of no 
more than ~20 Myr (see Dohm-Palmer & Skillman 2002). This implies that these stars 
should still lie very near their birth locations. This makes the distribution of supergiant 
candidates in the outer disk of IC 2233 noteworthy, as they show no obvious concentration 
along the midplane, but rather a wide spread of heights. This trend is particularly striking 
compared with the steep, roughly exponential vertical light profile seen near the central parts 
of the galaxy and implies an overall flaring of the young stellar disk. We also see evidence for 
flaring in both the ionized and neutral layers in IC 2233 (see §3.31 fc 1^^ . Interestingly, the 



^For comparison of our exponential disk fits with the isothermal disk fits of Seth et al. 2005, we assume 
hz ~ 0.549zo, where zq is the isothermal scale height. 
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scale heights of the gaseous and stellar components appear to be roughly comparable in the 
outer disk of IC 2233 — a factor that may be linked with the rather low global star-formation 
efficiency in this galaxy (see also Ferguson 1998; Ferguson et al. 1998; Matthews & Wood 
2003). 

3.3. The Ionized Gas in IC 2233 

Our continuum-subtracted Hq;+[Nii] image of IC 2233 is reproduced in Figure [2l Al- 
though our bandpass was broad enough to encompass the [N ll] lines at 6548A and 6583A in 
addition to Ha, the [Nil] emission in IC 2233 is known to be quite weak (Goad & Roberts 
1981) and therefore should contribute minimally to this image. 

The ionized gas distribution in IC 2233 has been studied previously by Rand (1996), 
who drew attention to plumes of emission extending to more than 1 diQ kpc from the plane, 
near the center of IC 2233 (see also Miller & Veilleux 2003). These features are also visible in 
our Figure [21 Rand (1996) suggested that this vertically-extended emission may be due to an 
outflow from a central starburst. However, comparison with the accompanying broadband 
image in Figure [2] shows that the region from which the ionized plumes originate corresponds 
(in projection) to a knotty, extremely blue [{B—R)o ~ 0.3] clump of starlight that is displaced 
from both the galaxy midplane and from the isophotal and kinematic centers of the galaxy. 
This region has an extent of ~ 380 diQ pc. Its structure, size, and color are consistent with 
a cluster of multiple OB associations (Pietrzyhski et al. 2005). The off-planar appearance 
of this "superassociation" and the lack of any other comparable features within the galaxy 
suggest the possibility that this could be a remnant from a dwarf intruder. This would be 
consistent with the very different emission line ratios found by Goad & Roberts (1981) for 
the central regions of IC 2233 compared with the southern end of the disk, although we flnd 
no kinematic evidence to support this hypothesis (see below; see also Garrido et al. 2004). 

As also noted by Rand (1996), IC 2233 does not have a widespread extraplanar layer of 
ionized gas, consistent with the low global star formation rate of this galaxy. Nonetheless, as 
seen in Figure 121 the ionized emission exhibits a great deal of complex structure, ranging from 
small-scale features (compact Hll regions, filaments, and shells) to diffuse planar emission, 
to bright, giant Hll complexes. 

The largest and most luminous H II complex visible in IC 2233 is situated on the southern 
end of the disk. This region has a projected extension of more than 2 dio kpc and contains 
~30% of the total Ha luminosity of IC 2233. The location of this complex well outside the 
central regions of the galaxy underscores that recent high-mass star formation in IC 2233 has 
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not been centrally concentrated. Giant, off-center Hll complexes are also seen in other LSB 
galaxies with similar masses to IC 2233, such as the nearby, face-on LSB galaxy NGC 4395 
(Wang et al. 1999). These complexes may be more prevalent in these low-mass disks with 
slowly rising rotation curves, since these systems lack the shearing forces that would tend to 
rapidly stretch and distort similar large-scale density enhancements in more massive spirals 
(Elmegreen & Hunter 2000). In addition, the weaker self-gravity and intrinsically thicker gas 
disks of small spirals may allow in general the formation of larger bound condensations (Kim 
& Ostriker 2006). Contrasting with these giant Hll complexes, some lines of sight through 
the IC 2233 disk show little Ha emission. This implies that the recent star formation in 
IC 2233 has been patchy and highly localized, perhaps occurring along loose spiral arms 
analogous to those visible in NGC 4395. 

Turning to smaller scales, several shell- or bubble-like features of various surface bright- 
nesses can be seen in the ionized ISM of IC 2233. One example is the circular, low surface 
brightness feature centered at aj2ooo.o = 08'^13""58.5^ 5j2ooo.o = +45°44'21'a (Figure [2]), 
which has a radius of ~ 2" (r ~ 100 dio pc) and a projected height of ~ 100 dio pc from 
the disk midplane. Its size and morphology are consistent with a wind-blown superbubble 
produced by an association of hot, young OB stars (Mac Low & McCray 1988). Supporting 
this interpretation, our B and R continuum images show a blue point source at the apparent 
center of the bubble. 

Another interesting aspect of the ionized gas morphology of IC 2233 is that the various 
H II regions and complexes in IC 2233 are not confined to a thin layer. This effect becomes 
more pronounced toward the outer parts of the galaxy, where a number of bright H ll regions 
or complexes appear to be displaced by ~ 0.5hz or more from the midplane. The brightest 
Hll complexes also seem to trace a rather wavy pattern about the midplane, although we 
are unable to identify a clear correspondence between these undulations and the more well- 
defined corrugated pattern seen in the neutral hydrogen disk ( §6.41) . Overall, the complex, 
"frothy" Ha morphology of IC 2233 forms an interesting contrast to the much flatter ionized 
gas disk of the superthin galaxy UGC 7321, which shows very little complex structure and 
no sign of giant (kiloparsec-scale) Hll complexes (cf. Figure 7 of Matthews et al. 1999). 

Although our Ha image is not flux-calibrated, we have estimated an Ha luminos- 
ity for IC 2233 by using our i?-band photometric calibration and the response curves for 
the R-hand (continuum) and narrow-band (Ha+[Nll]) filters. In this manner we derive 
Lua ~ 6.1 ± 1.8 X lO^^d^Q erg s~^. This value is corrected neither for internal extinction nor 
for the contribution from [Nil] emission in the narrowband filter (likely to be a ~10% effect; 
see Goad & Roberts 1981). Our estimate for Lhq is in good agreement with a value derived 
from the data of Rand (1996) (Lho ~ 5 x lO^^ci^Q erg s~^, with an uncertainty of ~ 10%; 



- 13 - 



R. Rand 2003, private communication). Using the relations given by Kennicutt (1998), the 
observed Ha luminosity translates to a global star formation rate of ~ O.OSMq yr"^. Thus 
the disk-averaged star formation rate in IC 2233 is quite low despite various signatures of 
ongoing, localized star formation (e.g., Hll complexes; a blue global color, OB associations, 
and candidate populations of supergiant stars). For comparison, we can also compute an 
estimate of the star formation rate based on the IRAS FIR emission. Using the 60 fim and 
100 /im fluxes from the NED database we find Lpm = (7.1±0.9) x W'^df^ Lq. Following Ken- 
nicutt (1998), this translates to a global star formation rate of 0.02 Mq yr^^, where we have 
assumed that the total IR luminosity is ~ 1.75Lfir. This estimate is more than a factor of 
two smaller than the Ha estimate, consistent with the suggestion of Bell (2003) that the FIR 
may systematically underestimate the total star formation rate of low-luminosity galaxies. 
In any case, both estimates represent the extreme low end of observed star formation rates 
for Sd spirals (Kewley et al. 2002). This is consistent with the low emission line intensity 
ratios measured in IC 2233 by Goad & Roberts (1981) and Miller & Veilleux (2003) and 
with the general classification of IC 2233 as an LSB spiral. 



4. VLA Observations 

We observed IC 2233 and NGC 2537 simultaneously in the H I 21-cm line using the VLA 
in its C (3 km) configuratiorfl on 2000 May 28 (hereafter Day 1) and 2000 May 29 (Day 2). 
The observations were performed over two 8-hour sessions, with a total of ~ 11.6 hours of 
data acquired on the target field. Approximately 97% of the data were obtained at elevations 
greater than 40°. 

Our observations were carried out using the 4-IF spectral line observing mode with a 
1.5 MHz bandpass per IF pair (left and right circular polarization) and on-line Hanning 
smoothing. This resulted in 63 channels of width 24.4 kHz (~5.2 km s~^) per IF pair. One 
IF pair was centered on the published redshifted velocity of IC 2233 (Vhei=565 km s~^) while 
the other IF pair was centered on that of NGC 2537 (\4ici=447 km s~^). The field center was 
chosen to be 5' north of the optical center of IC 2233 to prevent excessive attenuation at the 
position of NGC 2537, with only a ~ 7% loss at the position of IC2233. Five-minute observa- 
tions of J0713+438 were interspersed with 30-minute observations of the IC2233/NGC2537 
field in order to determine phase calibrations as well as bandpass corrections (see §4.1.ip . 
In addition, J0137+331 (3C48) and J1331+305 (3C286) were each observed for 10 minutes 



^Since 1998 the C configuration has been modified from its original form by placing the fifth antenna of 
the north arm at the center of the array in order to improve sensitivity to extended emission. 
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each day for use as flux calibrators. Details of our observations are summarized in Table 2. 

4.1. Data Reduction 

4.1.1. Calibration 

A vector average of the visibilities from channels 5-59 at each time stamp was used 
to identify interference and corrupted data. Because our observations were obtained in 
the daytime during solar maximum, our data contain significant solar interference at short 
spacings, and noticeably poorer phase stability compared with Hi observations of other 
targets obtained by us on adjacent dates, but after sunset (Paper I). The Sun was ~ 52° 
away from our target field during the observations and relatively "quiet" on our first day 
but far more active on the second day, which included a period of about one hour where 
only baselines between the farthest two antennas of each arm gave usable data. In addition, 
on the second day, a significant fraction of the data from three antennas (located at the 
East-2, East-8, and East- 16 stations, respectively) had to be discarded due to electronic 
problems. In spite of these losses, we have succeeded in making images of excellent quality, 
and comparison with past single-dish observations suggests that we have recovered the bulk 
of the Hi emission in both IC 2233 and NGC 2537 (see gEl&EH)- 

The flux density scale was determined using the primary calibrators 3C48 and 3C286 
and adopting fluxes of 15.90 Jy and 14.74 Jy, respectively, at 1417.64 MHz. These fluxes 
were scaled to our other frequency setting following the precepts of the VLA calibration 
manual. Using one of the two primary flux calibrators to solve for the flux of the other, 
and comparing this flux with its known (predicted) value, indicates that our internal flux 
calibration uncertainty is less than 1%. 

We determined antenna-based amplitude and phase calibrations using the observations 
of J0713+438 and an initial bandpass calibration from the average of the data obtained 
on 3C48 and 3C286. This initial bandpass calibration resulted in some curvature in the 
spectra taken through strong point sources in the IC 2233/NGC 2537 field. We therefore 
evaluated and applied further scan-based bandpass corrections by fitting 5th order Chebyshev 
polynomials to channels 5-59 of each (scan) averaged spectrum of J0713-I-438 and corrected 
the IC2233/NGC2537 data by linear interpolation (in time). This procedure led to flat 
bandpasses through the bright continuum sources in the field. The spectrum of the strongest 
source (4C-I-46.17) shows no significant deviations at the level of 0.13% (2.5 a). 

The IC 2233 field contains significant background continuum emission (totaling ~ 1.1 Jy 
observed, i.e. before correction for the attenuation of the primary beam) which is dominated 
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Table 2. Summary of IC 2233/NGC 2537 H I Observations 



Parameter Value 



Observing Set-Up 



Array configuration 
Baseline range 

Observing dates (2x8 hours) 
Pointing and phase center, a (J2000.0) 
Pointing and phase center, S (J2000.0) 
Total on-source observing time 
Flux calibrators 
Phase calibrator 
Number of IFs 

Channel width (after Hanning smoothing) 
Velocity separation of channels 
Number of channels per IF 

Heliocentric central velocity (optical definition) 
Primary beam (FWHM) 



C 

30 - 3385 m 

2000 May 28 (Day 1), 29 (Day 2) 
08 13 58.8 
+45 49 36.3 

696 minutes (683 min. at elevation > 40°) 
J0137+331 (3C48), J1331+305 (3C286) 
J0713+438 
4 [2 X (RR, LL)] 
24.4 kHz 
5.2 km s"-*- 

63 

565.0 km s^i (IF pair 1); 447.0 km s^^ (IF pair 2) 
~ 30.'6 



Deconvolved Image Characteristics 



Robustness parameter (7^)^ 

Synthesized beam FWHM'' 

Synthesized beam position angle 

Linear resolution of synthesized beam'' 

rms noise per channel (Icr) 

rms noise in total-intensity image (Icr) 

Limiting column density per channel (Icr) 

rms noise in Hi mass channel"^, Mjji (Icr)'' 



1 

~ 16'.'1 X 14'.'8 

70° 

~ 0.8 did kpc 
(0.37-0.44) mjy beam"' 
~ 5 mJy beam~^ km s~^ 
(8.8 - 10.5) xlO^* atoms cm-^ 
(4.5 - 5.3) X 10"' dlgMQ 



^We also used TZ = 0.7 (and a circular restoring beam with FWHM = 15'.'0) to make the continuum 
image (see i|4. 1.31 1 

^diQ is the distance to the target galaxies expressed in units of 10 Mpc (the distance to IC 2233, see 
Table 1). 
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by two bright point sources (4C+45.15 at ctjsooo.o = 08'*12"*42.0^ 5j2ooo.o = +45°36'51':3 

and 4C+46.17 at ajaooo.o = 08'^14™30.4^ 5j2ooo.o = +45°56'39"4). For these "4C" sources 
we observed flux densities of 531±5 mJy (best fit as two extended components with formal 
uncertainty below the 1% adopted here) and 1089±11 mJy (1% calibration uncertainty, with 
a formal error of 0.5 mJy). These measurements are in good agreement with values from 
the NRAO VIA Sky Survey (NVSS; 528 ± 16 mJy and 1106 ±33 mJy respectively; Condon 
et al. 1998). The strength of these two sources results in significant ghosts in the images 
of channels on both edges of each IF pair. These ghosts, which appear at the symmetric 
positions of these sources with respect to the phase center of the observations were first 
discussed by Bos (1984,1985) in his study of the Westerbork array. They are always present 
but rarely noticed at the VLA. We observe them at a strength of ~ 4.3% in channel 1 
and can trace them up to ~ channel 10 (although they would be difficult to detect beyond 
channel 4 if their spatial location were unknown). On the other end of the band, the ghosts 
are clearly seen in channels 61-63. The cumulative effect of such ghosts is responsible for 
the well-known increase in the noise of the channel images at the edges of the band in VLA 
spectral observations (Uson 2007). 

H I emission is observed in channels 13-63 of IFl and 1-45 of 1F2, leaving only a small 
number of line-free channels. Avoiding the ghost emission discussed above, we used chan- 
nels 5-12 of IFl and channels 46-59 of IF2 to improve our gain calibration through self- 
cahbration. We took the image obtained from channels 5-12 of IFl and day 1 as a reference 
to compute phase corrections in 10-minute intervals for the observations of both IFs and 
both days (using channels 5-12 for IFl and channels 46-59 for IF2). After applying the 
corresponding phase corrections, we determined amplitude and phase corrections with a 30- 
minute averaging time (one set per scan) restricting the baselines used to those larger than 
200A (the shorter baselines contained significant amounts of corrupted data; sec above). The 
average of the amplitudes was constrained to be constant in order to preserve the overall 
amplitude calibration. 

The frequency-averaged residual visibilities (after subtraction of the theoretical contri- 
bution of the detected continuum emission) for each day and each IF were used to find 
low-level contamination of the data. We discarded the corresponding data and performed 
one additional iteration of phase-only self-calibration on 10-minute intervals which led to 
improved pseudo-continuum images for each day and IF. 
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4.I.2. Imaging of the Line Emission 

The strong continuum present in our field precluded a fast subtraction of the continuum 
emission by fitting the real and imaginary components of each visibility (Cornwell et al. 
1992). Instead, we have used an alternate, two-step continuum subtraction procedure. First, 
we computed a continuum model for each IF pair on each day by imaging the available line- 
free channels. We used channels 6-12 for IFl and channels 46-59 for IF2. The data were 
gridded channel by channel to avoid bandwidth smearing. We subtracted these models 
(channel by channel) from each respective database and concatenated the residual data for 
each IF pair for each of the two days. Next, we obtained channel images using a 3" cell 
size and a robustness factor TZ = 1. The resulting synthesized beam was 16''0xl4''8 at 
a position angle of —70°. Deconvolution to a level of 0.40 mJy beam~^ was followed by 
restoration of isolated components if their absolute value inside a 6-pixel radius was less 
than 1 mJy beam~^. The rms noise in the resulting image cube ranged from 0.40-0.46 mJy 
beam~^ channel"^ (for channels 5-59). We found that this procedure removed the bulk of the 
continuum emission, although some low-level sidelobes (of order ±1 mJy) from the brightest 
continuum sources in the field remained. 

To improve the continuum subtraction further, we combined the data for the two IFs 
(which overlap in frequency) to make a "global cube" containing 86 channels and a 2.1 MHz 
bandwidth. The corresponding channels from both IFs were averaged with equal weight 
after discarding four edge channels of each IF (60-63 of IFl and 1-4 of IF2). The frequency 
misalignment between the two IFs was approximately one-fifth of a channel, so this resulted 
in a small amount of frequency smearing (~1.1 km s^^); however, this is inconsequential 
for the analysis of the spectral images that follows. We blanked regions that contained 
line emission in the global cube (leaving only the residual continuum) and fit a first-order 
polynomial (in frequency) to each (spatial) pixel in the blanked cube (see Cornwell et al. 
1992). Finally, we used the average level and slope derived from this fit to remove the 
residual continuum from each pixel/channel of the global cube containing the line emission. 

The rms noise in the resulting continuum-subtracted global cube was reduced to 0.37- 
0.44 mJy beam""^ channel"^, and the noise showed a Gaussian distribution exclusive of the 
line emission. Noise in the channels where the data from the two IF pairs were averaged 
have ~10% lower noise. Spectra through the locations of the brightest continuum sources 
appear extremely flat, with an rms noise of ~0.42 mJy and mean flux density levels across 
the bandpass that are indistinguishable from zero. We have used this resulting continuum- 
subtracted global image cube for the analysis of the line emission in IC 2233 and NGC 2537 
described hereafter. 
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4.1.3. Imaging of the Continuum Emission 

Although the continuum images used to process the spectral line data were adequate for 
the purpose described in the previous section, they contained a number of systematic errors 
that became more pronounced when we averaged the four images (2 IFs, 2 days) to obtain 
a higher sensitivity continuum image for the purpose of studying the emission of IC 2233 
and NGC 2537. As discussed in §4.1.11 the field contains two "4C" sources that are located 
at positions where the response of the primary beam is at ~ 80% and ~ 35%, respectively, 
of its maximum value. In addition, the off-axis geometry of the Cassegrain optics of the 
VLA antennas induces a beam-squint which separates the right and left circular polarization 
beams by 104" on the sky (Weinreb et al. 1977). The response at the position of these 
sources is different for both polarizations and it is modulated oppositely with parallactic 
angle. Therefore, the visibilities do not obey a simple convolution equation and the image 
contains a significant and serious error pattern in the form of concentric rings around the 
two main continuum sources. This pattern extends over the full field-of-view and affects the 
flux density at the locations of IC 2233 and NGC 2537. 

Although the beam-squint imposes a significant modulation on the continuum emission 
in the IC 2233/NGC 2537 field, the effect is negligible provided that the data are restricted 
to those visibilities for which both polarizations are available. In addition, it is important to 
ensure that amplitude self-calibration does not "follow the squint" on the strongest source 
in the field. Thus, we have made an image of the continuum emission, restricting the 
data to those visibilities for which both polarizations are available and pre-averaging both 
polarizations prior to determining corrections to the amplitudes through self-calibration. We 
have restricted the baselines to those longer than 1600A to minimize low-level cross-talk and 
solar contamination. This results in a 12% drop in the number of visibilities compared to the 
previous estimate of the continuum ( §4.1.11) . Furthermore, we have found that the ghosts 
discussed in §4.1.11 are still present in channels 5 and 6 of IFl and in channel 59 of IF2. 
Therefore we made a new image using channels 7-12 of IFl as well as channels 46-58 of 
IF2. We used "3-D" gridding and a robustness factor 71 = 0.7, which yields the optimal 
compromise between noise and resolution for VLA imaging observations obtained with full 
u-v coverage (Briggs 1995). The resulting image has an rms noise of 0.12 mJy beam~^ and 
is free of systematic errors, except at the positions of the two brightest "4C" sources, where 
it is still somewhat dynamic-range limited. This image, which corresponds to an average 
frequency of 1418.3 MHz, is shown in Figure [31 
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5. The 21-cm Continuum Emission in the IC 2233 Field 

As previously described, the IC 2233/NGC 2537 field contains significant continuum 
emission from background sources. The observed fiux density (uncorrected by the attenua- 
tion of the primary beam) in the central 1° x 1° is 1.129 Jy, with a total of 1.145 Jy when two 
flanking fields located outside the main beam are included. The fiux density is dominated 
by the two aforementioned "4C" sources (see §4.1.ip . 

For NGC 2537 we detect continuum emission totaling 10.3±0.6 mJy — in good agreement 
with the NVSS (10.5±1.7). The peak emission of 2.4 ±0.2 mJy lies at ajaooo = 08'^13™13.1^ 
^j2ooo = +45°59'40", consistent with the location of the most prominent star forming region 
in this galaxy. 

We do not detect any significant continuum emission from IC 2233. Within an area 
defined by the stellar disk of the galaxy we measure an integrated fiux of 1.7 ± 0.6 mJy, 
after correction for the attenuation of the primary beam. However, given the dynamic range 
limitations and residual sidelobes discussed in the previous section, this should be considered 
non-significant and converted to a (2cr) upper limit of <3 mJy. The weak or absent radio 
continuum emission from IC 2233 is consistent with its LSB nature and with the low global 
star formation rates derived from other tracers, including Ha and FIR emission (see §3.3p . 
Using the FIR luminosity of IC 2233 from Table 1, the radio-FIR correlation of Condon 
(1992) predicts a 1.4 GHz continuum fiux of ~ 3 mJy, consistent with our upper limit. 

6. Properties of the Neutral Hydrogen in IC 2233 

6.1. Channel Images 

In Figure H] we show the individual channel images for all velocities where H I line emis- 
sion was detected in IC 2233. Several key properties of the Hi distribution and kinematics 
of IC 2233 can be gleaned from these images. 

There are significant differences in the distribution and morphology of the H I emission 
on the approaching versus the receding side of IC 2233. For example, on the receding 
(northern) side the emission distribution in the individual channels appears more radially 
elongated; all but the outermost few channels show emission near a; = (where we take x to 
lie along the major axis of the galaxy, with positive x toward the north), and the channels in 
the velocity range 627.1-637.4 km s~^ show emission extending along the full extent of the 
receding side of the disk (x ~ — 200"). Within these latter channels, the peak brightness 
occurs ~ 80" from the galaxy center. On the approaching (southern) side, the corresponding 
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channels show a more compact distribution of emission, with brightness peaking ~ 120" from 
the galaxy center. These differences imply an asymmetry in the H I distribution and velocity 
field on the two sides of the galaxy (see also below). 

The channel images in Figure H] also provide insight into the vertical structure of the neu- 
tral hydrogen in IC 2233. On both the approaching and receding sides of the galaxy — most 
notably in intermediate velocity channels — the isointensity contours broaden with increas- 
ing projected distance from the galaxy center, implying an overall flaring of the gas layer 
with increasing galactocentric radius (see also §6.2p . Three-dimensional modeling would be 
necessary to characterize the flaring in detail (see Oiling 1996; Matthews & Wood 2003). 
However, because IC 2233 is nearly exactly edge-on, measurements of the vertical extent of 
the emission in select channel images can provide estimates of how much the gas layer thick- 
ness changes as a function of radius. Channels near the velocity extrema contain emission 
predominantly from the outer regions of the disk. Therefore, the thickness of the emission 
along the z direction in these channels provides an estimate of the maximum thickness of 
the gas layer. At the same time (owing to the curvature of the isovelocity contours), emis- 
sion observed near a; = in channels of intermediate velocity should arise from locations 
close to the true galaxy center and be largely uncontaminated by outer disk gas along the 
line-of-sight. The vertical thickness of this emission therefore should give an estimate of the 
H I scale height in the inner galaxy. 

To gauge the thickness of the H I layer in the outer disk of IC 2233, we have used Gaussian 
fits along the z direction to the emission in channels corresponding \,oV = 461.6 km s~^ and 
V = 647.8 km s~^. These channels provide better signal-to- noise than the very outermost 
channels, but still sample primarily outer disk emission. The fitted intensity profiles were 
extracted through the location of the peak intensity in each of these channels. The resulting 
deconvolved FWHM thickness in both cases is ~ 20.3" (~ 1 diQ kpc). To estimate the 
thickness of the inner gas disk, we limit our measurement to the receding side of the galaxy. 
We fit the emission along x ~ in the channel corresponding to \^=606.4 km s~^ and find 
a deconvolved FWHM of 10"4 (~ 500 dio pc). Similar results were obtained in adjacent 
channels. We conclude that the H I scale height of IC 2233 increases by at least a factor of 
two from the center to the outskirts of the galaxy. 

In addition to the trends described above, several channels show plumes of emission 
stretching to heights of > 50" above and below the plane that cannot readily be attributed 
to warping or flaring of the disk. This phenomenon is particularly prevalent on the receding 
side of the galaxy. For example, in the channel images corresponding to heliocentric radial 
velocities between ~ 586 — 658 km s~^, we see a faint envelope of "extra" emission lying 
between roughly x = 100" — 150", on both the ±z sides of the disk. The high- 2; material shifts 
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little in position from channel to channel; therefore toward the highest velocity channels, this 
gas lies at smaller projected radii than most of the emission, while in channels corresponding 
toV = 586 — 617 km s~^, this material is seen "trailing" behind the bulk of the Hi. Although 
there is relatively little H I mass associated with these extraplanar features, the large spread 
of velocities observed near a single spatial location implies a large velocity dispersion for this 
material. The characteristics of this emission are therefore consistent with those expected for 
a rotationally lagging Hi "halo" (see Figure [12] below), analogous to those now identified in 
several other spiral galaxies (e.g., Swaters et al. 1997; Fraternali et al. 2001), including the 
LSB spiral UGC 7321 (Matthews & Wood 2003). The possible presence of an Hi halo around 
IC 2233 is of particular interest as a second example of this phenomenon in a galaxy with a 
low global star formation rate and minimal extraplanar ionized gas component. However, a 
more rigorous characterization of this extraplanar H l emission will require three-dimensional 
kinematic modeling and is beyond the scope of this paper. Note that the estimates of disk 
flaring computed above should be largely unaffected by a rotationally lagging halo, since the 
halo gas is expected to manifest itself as non-Gaussian wings to the vertical intensity profiles 
rather than a global broadening (see Matthews & Wood 2003). 



6.1.1. The Kinematic Center of IC 2233 

The systemic velocity and kinematic center of spiral galaxies often are defined using 
either the centroid of the global H I profile or the location that minimizes the asymmetry in 
the rotation curve. However, given the frequency of lopsidedness and other perturbations 
in galaxies (see §9.11) . such choices are not always physically justified. Here we propose an 
alternate method of defining the systemic velocity and kinematic center of an edge-on galaxy. 
Our method is most appropriate for late-type spirals, where Hi typically has a high filling 
factor. 

For a differentially rotating disk whose velocity field is spatially and spectrally resolved, 
the gas emission will be expected to have its minimum observed spatial extent along the 
direction of the disk major axis (Axmin) within the channel corresponding to the velocity 
range V^ys±0.5AV^, where V^js is the true systemic velocity and is the spectral resolution. 
For all other observed velocities, V^, emission observed within V^iO.SAV^ will span a larger x- 
extent (Axj > Axmin) owing to the projection along the x direction of the curved isovelocity 
loci (see e.g.. Figure 6 of Sancisi & Allen 1979). In the absence of severe optical depth effects, 
the emission is also likely to exhibit its peak brightness in the channel encompassing the true 
systemic velocity owing to the combination of its minimum projected spatial extent and the 
maximum path length through the disk. 
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Using elliptical Gaussian fits, we derived a position angle and FWHM major and minor 
axis diameters for the Hi emission distribution in each of the channel images of IC 2233 
(corrected for the broadening due to the synthesized beam). A simple parabolic fit to the 
values of the major-axis diameters (Axj) as a function of channel yielded a minimum at 
channel 34.25 ± 0.1 with ~ 5 with 5 degrees of freedom (Figure [5]), corresponding to a 
heliocentric velocity of 553.4 ±0.5 km s~^. The orientation of the major axis in channel 34 
corresponds to a PA= 172.0° ± 0.9° , well aligned with the disk of IC 2233. A similar series 
of fits to the peak H I brightness as a function of channel yielded a maximum corresponding 
to channel 34.75 ± 0.2. However, self-absorption as well as local changes in the gas phase 
can affect the column density of H I . In addition, the peak brightnesses in channels 34 and 
35 are essentially identical. Therefore, we adopt a value of V^ys = 553.4 ±1.0 km s^^ (2cr) 
for the systemic velocity of IC 2233. We derived the location of the kinematic center of 
IC 2233 by linear interpolation to "channel" 34.25 using the centroids of the Hi emission in 
channels 34 and 35, which yields: aj2ooo.o = 08'^13"^58.9^ (5j2ooo.o = ±45°44'27':0 with formal 
uncertainties of 0.1** in right ascension and 1''0 in declination (2(t). The location of the Hi 
kinematic center coincides within the errors with the location of the minor axis dust feature 
seen in our optical images of IC 2233 (see §3.2j Figure [2]). 



6.2. The Total Hi Content 

To derive a global Hi profile for IC 2233 and measure its integrated Hi flux, we have 
applied a "percolation" procedure to the individual channel images containing galactic line 
emission. For each channel we defined an irregular blotch containing line emission above 
1 mJy beam^^ (~ 2.5cr) and successively expanded this region for n iterations using bands 
of 2 pixels, until the flux within the blotch converged. Each step increases the noise in the 
total flux, as it is proportional to the square-root of the number of synthesized-beam areas 
inside the blotch. Convergence in each channel was assumed if an expansion band lead to 
no significant increase in the total flux (i.e., the SNR of the total flux did not increase with 
an additional step). To avoid biases, the algorithm begins by finding pixels with absolute 
value above the 1 mJy beam~^ threshold, with spurious noise spikes (positive and negative) 
being discarded. We corrected the channel images for primary beam attenuation prior to 
determining the flux in the regions localized with the percolation algorithm, and estimated 
an error based on the rms noise in the channel and the number of beam areas in the blotch. 
These uncertainties are <4 mJy channel"^. The resulting global Hi profile is shown in 
Figure [61 

The global H I profile of IC 2233 exhibits a classic double-horn shape, but with a clear 
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lopsidedness. As illustrated in the inset of Figure O the slopes of the two edges of the profile 
are also different: the side having the brighter horn shows a steeper slope than the side with 
the weaker and narrower horn. We comment on the significance of this further in §9.11 

From our global H I profile we measure velocity widths at 20% and 50% of peak maximum 
of 195.0 km s~^ and 174.1 km s~^, respectively. These values are in good agreement with 
previous measurements in the literature (e.g., Haynes et al. 1999). The centroid of our global 
H I profile as defined by the 20% peak maximum value is 557.1 km s~^. This is slightly higher 
than the systemic velocity of 553.4±1.0 km s~^ defined by the kinematic center of the galaxy 
(DID- 

The integrated Hi flux that we derive for IC 2233 is 47.2 ± 0.6 Jy km s'^. The error 
includes a statistical contribution from each channel (between 1 mJy and 4 mJy) and a 
calibration uncertainty of ~1%. Our integrated flux is consistent with the mean of various 
published single-dish values including: 53.1±4.6 Jy km s~^ (Fisher & Tully 1981); 43.1±2 Jy 
km s-i (Tifft 1990); and 47.55±0.5±2.0 Jy km s-^ (Haynes et al. 1999). Published values of 
the integrated H I flux from aperture synthesis measurements include 44.9±2.0 Jy km s~^ (Stil 
& Israel 2002a) and 52.3±7.8 Jy km s^^ (Swaters et al. 2002). Assuming that the H I emission 
in IC 2233 is optically thin, we derive a total Hi mass of Mhi=(1.11±0.01) x 10^ (IIq Mq. 



6.3. The Hi Intensity Distribution 

We derived the total Hi intensity distribution (zeroth moment map) for IC 2233 by 
summing the emission identified in each channel by the percolation algorithm described in 
the previous section. The resulting total H I intensity map is shown in Figure [71 overplotted 
as contours on our Hq;+[N ii] image. Our new H I map of IC 2233 reveals a variety of features 
not seen in previous H I images of this galaxy. 

Along the midplane of IC 2233, the H I distribution of IC 2233 appears clumpy in spite 
of the line-of-sight averaging effects of our edge-on viewing angle. At the resolution of our 
data, these clumps have peak column densities of ~ (6 — 8) x 10^^ atoms cm~^ — roughly 
30% higher than the surrounding gas. The clumps are not distributed symmetrically about 
the kinematic center of the galaxy, nor do they lie in a flat plane (see also §6.4p . As seen 
in Figure [71 the locations of all but the northernmost clump correspond with some of the 
brightest H ll complexes in the galaxy. Viewed face-on, the H I disk of IC 2233 may appear 
similar to that of the dwarf LSB spiral NGC 4395, where the H I shows a lopsided distribution 
with little central concentration and a number of denser, sub-kiloparsec-scale clumps are 
visible scattered throughout the disk (see Swaters et al. 1999). Most of these latter clumps 
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correspond with the locations of the brightest star-forming regions in NGC 4395. 

Further from the midplane of IC 2233, the H I isophotes become increasingly complex 
in shape. Various ripples and protrusions form a network of material extending to \z\ < 60" 
along the full radial extent of the galaxy. The bulk of these features cannot be attributed to 
noise, as most are visible in multiple successive contours. Figure [7] illustrates that the full 
vertical extent of the H I disk of IC 2233 is significantly larger than that of the ionized gas 
disk, even after accounting for resolution effects. While part of this effect comes from flaring 
and warping of the Hi layer, it seems likely that at least some portion of this vertically- 
extended emission comprises a distinct Hi "halo" (see also §6.11 and §7.21) . 

Along the radial direction, we measure the Hi extent of IC 2233 to be 6'66±0'04 at a 
limiting column density of 10^° atoms cm~^. This is ~1.3 times the optical diameter (defined 
by the 25.5 mag arcsec^ i?-band isophote). Perpendicular to the disk, the deconvolved 
FWHM thickness of the H I layer through the location of the kinematic center is 20''0±0''2 
(~ 970 dio pc). Similar fits to the emission in the individual channel images ( §6.ip show 
that the Hi layer in the inner galaxy is intrinsically thinner, then flares at larger radii, 
especially on the northern side where the (projected) thickness of the gas layer reaches 
FWHM ~ 30". Despite projection effects, this flaring is also evident in Figure [71 where 
we find that the global, deconvolved FWHM of the Hi layer increases systematically with 
increasing galactocentric radius. 

It is interesting to contrast the vertical structure of the Hi disk of IC 2233 with that of 
the superthin galaxy UGC 7321 studied in Paper I. We measured the global FWHM of the 
H I disk of UGC 7321 to be ~ 810 rfio pc — roughly 15% smaller than we measure for IC 2233. 
In addition, the distance- independent Hi axial ratio of UGC 7321 is significantly larger than 
that of IC 2233: {a/b)^! ~29 versus ~18, respectively. These differences cannot be attributed 
to projection effects given the similar orientation of the two galaxies, and imply that the 
mean thickness of the Hi disk of IC 2233 is intrinsically larger than that of UGC 7321, 
consistent with its thicker stellar and ionized gas disks. While part of this difference could 
result from dynamical heating due to a recent minor merger (e.g., Reshetnikov & Combes 
1997) we find no direct evidence of such an event (see §9.ip . and it is quite hkely that at least 
part of the difference in thickness stems from the lower total mass of IC 2233 and thus the 
lower self-gravity of its disk. As the H I velocity dispersion of galaxy disks seems to be nearly 
constant at ~6-9 km s~^ in disk galaxies of a wide variety of types and morphologies (van 
der Kruit & Shostak 1984; Dickey 1996), it follows that the gas scale height of lower- mass 
disks should be proportionately larger (see also Brinks et al. 2002). This in turn may have 
important implications for the formation of molecular clouds in low-mass galaxies and their 
overall regulation of star formation (Elmegreen & Parravano 1994; Ferguson 1998; Dalcanton 
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et al. 2004; Matthews et al. 2005). 

Knowledge of the gas surface density of a galaxy is important to characterizing its star- 
formation efficiency. To estimate the deprojected Hi surface density in the disk of IC 2233, 
we have used an Abel inversion technique (e.g., Binney & Tremaine 1987). We extracted 
a 3"- wide slice along the major axis of the galaxy, sampled at 16" increments, and then 
computed deprojected intensity distributions for the two sides of the galaxy independently, 
using a fortran program adapted from Fleurier & Chapelle (1974). The results are presented 
in Figure [HI We assume that the H I is optically thin and cylindrically symmetric on each 
of the two sides of the galaxy. Because the observed intensity profile fiuctuates rapidly near 
the galaxy center, the deprojected surface density at small radii (r <30" or <1.5 dio kpc) is 
not well constrained. To minimize the impact of this effect, we have smoothed the measured 
intensity profile by a factor of three before each fit. Despite these small-scale uncertainties, 
the resulting profiles for the two sides of the galaxy should give a reasonable estimate for 
the mean H I surface density as a function of radius. Indeed, the integrated H I mass derived 
from the deprojected profiles shown in Figure [8] is 9.9 x 10* Mq — in excellent agreement with 
our observed value (see Table 3). Modest differences are apparent on the two sides of the 
disk. At all radii, the mean H I surface density of IC 2233 is considerably below the mean 
value of ~ IOMq pc~^ found by Cayatte et al. (1994) for high surface brightness Sd galaxies. 
However, values of Shi within the disk of IC 2233 are comparable to those seen in other LSB 
spirals (van der Hulst et al. 1993), including the superthin LSB spiral UGC 7321 (Paper I). 



6.4. Warping and Undulations in the IC 2233 Disk 

A careful examination of Figure [7] reveals that neither the brightest H I clumps nor the 
underlying, smoother gas component in IC 2233 lie along a flat plane. To illustrate this 
more clearly, we have plotted in Figure [9] the vertical displacement in the centroid of the 
H I layer as a function of radius. These displacements were derived from Gaussian flts to a 
sequence of slices perpendicular to the mean galactic plane. The resulting plot highlights 
two important features of the IC 2233 HI disk. 

At the edge of the stellar disk of IC 2233 (denoted by arrows on Figure [9]), we see that 
the Hi layer bends away from the plane, forming an "integral sign" warp. The amplitude 
of the warp is asymmetric, with an amplitude ~50% higher on one side of the disk than on 
the other {Az ~ +3''5 versus Az ^ —5".5). On the higher-amplitude side, the gas layer is 
observed to twist back toward the midplane at its outermost extent. These various traits are 
common features of galactic warps (Burton 1988; Garcia-Ruiz et al. 2002; Sanchez-Saavedra 
et al. 2003; Paper I). Recently Saha & Jog (2006) have suggested that amplitude asymmetries 



- 26 - 



Table 3. 21-cm Radio Properties of IC 2233 



Parameter 


Value 


Measured Quantities 


a (J2000.0)'' 


08 13 58.9 (±0.1) 


5 (J2000.0)^ 


+45 44 27.0 (±1.0) 


Peak H I column density'' 8 X 10^^ atoms cm ^ 




6.'66±0.'04 


eHi.b"^ (FWHM) 


22'.'6±0'.'2 


(a/6)Hi'= 


18 


/Fnidi. 


47.2 ±0.5 Jy km s"! 


W20* 


195.0 km s-i 


1^50* 


174.1 km s-i 




553.4±1.0 km s"! 




557.1 km s^^ 


Vmax 


~85 km s-i 


Fcont (21-cm) 


< 3 mjy (2(t) 


Derived Quantities 


Mm' 


(1.11 ±0.01) X 10^ dfo M© 




1.55x10^0 du) Mq 


Mhi/Ls 


0.62 Mq/Lq 


Mdyn/Mni 


14 d-i 


-DHl/i?25.5 


1.3 



''Kinematic center (see i|6.1.1|l 
''At 16" resolution. 

'^Measured at a column density of lO^'' atoms cm~^. 

'^Projected thickness of Hi layer along minor axis, mea- 
sured at a column density of lO'^" atoms cm~^ and cor- 
rected for the resolution of the synthesized beam. 

°Axial ratio of H I disk measured at a column density of 
10^" atoms cm~^ and corrected for the resolution of the 
synthesized beam. 

^Wp is the H I profile width measured at P% of the mean 
peak value of the two profile horns. 

^From kinematic method (heliocentric frame; optical def- 
inition). 

''From global Hi profile (heliocentric frame; optical def- 
inition). 

'Assuming the H I is optically thin and dio is the distance 
to IC 2233 expressed in units of 10 Mpc. 

jProm Mdyn = 2.326 X l0^rV'^{r), where we have taken 
r = 9.2 dio kpc (190") and V{r) =85 km s"!. 
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Table 4. IC 2233 Rotation Curve 



Major axis distance (arcsec) 


Circular Velocity (km s ^) 


Uncertainty (km s ^) 


-189.0 


472.8 


4.9 


-171.0 


473.9 


2.8 


-162.0 


470.5 


1.5 


-153.0 


472.9 


0.8 


-144.0 


473.0 


0.7 


-135.0 


475.7 


0.4 


-126.0 


467.0 


0.3 


-117.0 


469.2 


0.2 


-108.0 


475.9 


0.1 


-99.0 


484.0 


0.1 


-90.0 


488.8 


0.1 


-81.0 


491.5 


0.1 


-72.0 


493.2 


0.1 


-63.0 


498.0 


0.1 


-64.0 


504.0 


0.1 


-45.0 


508.7 


0.1 


-36.0 


514.5 


0.1 


-27.0 


523.4 


0.1 


-18.0 


531.7 


0.1 


-9.0 


542.6 


0.1 


0.0 


551.5 


0.1 


9.0 


573.6 


0.1 


18.0 


589.6 


0.1 


27.0 


597.3 


0.1 


36.0 


599.4 


0.1 


45.0 


603.9 


0.1 


54.0 


608.3 


0.1 


63.0 


617.0 


0.1 


72.0 


618.2 


0.1 


81.0 


620.3 


0.1 


90.0 


623.9 


0.1 


99.0 


627.8 


0.1 


108.0 


632.7 


0.1 


117.0 


635.7 


0.1 


126.0 


639.0 


0.1 


135.0 


645.8 


0.1 


144.0 


641.3 


0.1 


153.0 


641.8 


0.1 


162.0 


642.7 


0.2 


171.0 


641.6 


0.3 


180.0 


642.8 


0.4 


189.0 


640.0 


0.6 


198.0 


641.4 


1.1 


207.0 


641.2 


3.0 
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in warps are a natural consequence of dynamical wave interference between m=0 and m=l 
bending modes of a disk. The process that might excite such modes is still unclear, but one 
possibility is accretion of gas — either from a minor merger or from the intergalactic medium. 
We return to this issue in §9.1[ 

In addition to the warp, Figure [9] reveals another particularly intriguing feature of 
the neutral gas disk of IC 2233: a pattern of positive and negative vertical displacements 
spanning the full extent of the stellar disk. These undulations appear remarkably regular, 
with a wavelength of ~150" (7 diQ kpc). Their amplitude increases with distance from the 
center of the galaxy, reaching a maximum of \Az\ ~ 3" (~ 150 diQ pc). 

In the Milky Way, systematic deviations of the H I layer from the mean principal plane 
along the z direction (so-called "corrugations") were first noted by Gum et al. (1960). 
Subsequent studies found that corrugations are present along both radial and azimuthal 
directions over a significant fraction of the disk of the Galaxy (e.g., Quiroga 1974,1977; 
Spicker & Feitzinger 1986). The corrugations are reflected not only in the distribution of 
neutral hydrogen, but in a variety of tracers, including molecular clouds and Hll regions 
(Lockman 1977; Alfaro & Efremov 1996 and references therein). As measured in Hi, scales 
of the corrugations range from ~50-350 pc in amplitude and from one to several kpc in 
wavelength (Spicker & Feitzinger 1986), comparable in scale to the undulations seen in 
IC 2233. Although it has been suggested that corrugations should be a common phenomenon 
in disk galaxies, similar phenomena have so far been observed (in optical light) in only a few 
galaxies (Florido et al. 1991,1992; Alfaro et al. 2001). To our knowledge, corrugations have 
not been reported previously in the Hi disk of an external galaxy. 

Inspection of Figure [2] reveals that both the Ha emission and the optical continuum 
of IC 2233 also show deviations from a purely planar geometry, although these undulations 
are not as regular as those seen in Hi. Analysis of the vertical displacements via optical 
tracers is complicated by the patchy distribution of bright H ll regions and by the presence 
of a superimposed foreground star in the northern part of the galaxy. In the case of the 
i?-band data, we find a displacement of the intensity-weighted centroid from the midplane 
of Az =+l''8 at a; = —140"; the displacement systematically decreases with increasing x, 
reaching a value of Az = —2".7 at a; ~ —80", before increasing again monotonically to a value 
of Az ^0 near x = 0. The changes in Az with x thus have the same sense as in the H I data, 
but with different amplitude and phase. On the northern side of the disk, analysis is hindered 
by the bright foreground star, and no systematic trend in the displacement of the centroid as 
a function of radius is evident. All measured centroids between x = 50" — 125" have positive 
Az, with the exception of one point near x ~ 100" with Az ^0. The midplane deviations 
measured from the Ha data show no obvious systematic trend with radius, and using an 
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autocorrelation analysis we were unable to identify a periodicity in the Ha undulations on 
spatial scales > 10". These differences between the vertical displacements measured from the 
H I and optical data may result from intrinsic differences between the vertical displacements 
of the stars and the gaseous components, but projection effects and differences in the radial 
extent and filling factor of the different tracers are also likely to contribute. Disentangling 
these effects in an edge-on galaxy is further comphcated by the fact that disk corrugations 
may occur along the plane as well as orthogonally to it (see Spicker & Feitzinger 1986). 

The origin of the corrugations in the disk of IC 2233 is a puzzle. We will present a more 
extensive investigation of this issue based on additional multiwavelength observations in a 
future paper (L. D. Matthews & J. M. Uson, in preparation). Here we offer a brief discussion 
of this question based on results from the present data. 

In the Milky Way and in the face-on spiral galaxy NGC 5427, clear relationships have 
been found between corrugations and existing spiral arms, suggesting that the two phenom- 
ena are linked (Quiroga 1977; Spicker & Feitzinger 1986; Alfaro et al. 2001; see also Nelson 
1976; Martos & Cox 1998). While the edge-on orientation of IC 2233 makes it difficult to 
assess what type of spiral arm structure it might possess, the patchy distribution of Ha 
emission in this galaxy is suggestive of H ll regions organized along loose spiral arms. This 
is supported by the observed correlation between the locations of the Hi density enhance- 
ments along the midplane and the brightest Hll complexes (Figure [7]). The Ha rotation 
curve derived by Goad & Roberts (1981; see Figure \T3\ discussed below) also shows a series 
of amplitude "wiggles" similar to those linked with the spiral arms and velocity corrugations 
in face-on spirals (Alfaro et al. 2001 and references therein). 

Using numerical simulations, Edelsohn & Elmegreen (1997) showed that in a Milky Way- 
like galaxy, tidal perturbations from a low-mass companion will induce kpc-scale vertical 
oscillations (i.e., vertical displacements coupled with velocity undulations) with a spiral- 
like shape. In addition, simulations by Sellwood et al. (1998) have shown that vertical 
"buckling" — qualitatively similar in appearance to the undulations we see in IC 2233 — can 
be induced by the passage of a low-mass satellite through the disk. However, a problem 
with invoking either of these scenarios is that our optical and H I studies of IC 2233 have not 
uncovered any direct evidence of an ongoing tidal interaction or minor merger (see also ^ 
§9.ip . and Sellwood et al. found that if the satellite has been completely disrupted, bending 
waves will not be excited. While we noted in §3.21 that the unusual OB "superassociation" 
to the north of the kinematic center could be an intruder remnant, supporting kinematic 
evidence is lacking. Moreover, it is predicted that disk undulations would not persist once 
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the orbit of the satelhte has decayed into the disk (Sellwood et al. 1998)0 

If the midplane undulations in IC 2233 are not caused by either a recent intruder or 
hnked with (pre-existing) spiral arms, another possibility is that they reflect spontaneous 
bending instabilities in the disk (e.g., Griv & Chiueh 1998). This might offer a natural 
explanation for the increase in amplitude of the undulations with radius, as such an increase 
could be a consequence of the change in the group velocity of this wave as it reaches regions 
with larger gas scale-height and lower gas density in the outer disk (Spicker & Feitzinger 
1986; Hofner & Sparke 1994). If bending waves are the cause of the corrugations then similar 
features should be observable in the disks of other isolated galaxies. 

One final aspect of the corrugated structure of IC 2233 that is worth noting is that, 
regardless of its origin, this type of coherent, large-scale pattern can only exist if the gas disk 
of this LSB spiral is largely self-gravitating (e.g., Bosma 2002; Revaz & Pfenniger 2004). 
This adds to other recent evidence for "heavy" disks in some LSB disk galaxies based on the 
analysis of their spiral structure (Fuchs 2002,2003; Masset & Bureau 2003). 



7. The H I Kinematics of IC 2233 

7.1. H I Velocity Field 

The H I velocity field (first moment map) of IC 2233 is shown in Figure Near the 
midplane this diagram shows fairly straight and parallel isovels, as expected from the high 
inclination of the galaxy and its inner-disk kinematics, which resemble a solid-body rotator 
(see also §7.21 below). However, away from the midplane, many of the isovels curve to form 
S- or L-shaped lines. Such pronounced twisting of the velocity field is rarely observed in 
edge-on galaxies and could be a signature of warping or of gas on non-circular orbits (e.g., 
as a result of a lopsided or triaxial potential or large-scale streaming motions through spiral 
arms). Because IC 2233 is viewed so close to edge-on, the curved portions of the isovels trace 
primarily faint, extraplanar material, hence warping seems to be the most likely explanation, 
although the presence of some gas on non-circular orbits and/or a lopsidedness in the overall 
potential cannot be excluded (see also §9.11) . In addition, the kinematically anomalous gas 
discussed in §6.11 and §7. 21 may also contribute to the observed complexity of the extraplanar 
velocity structure. 



^We note that the simulations of Sellwood et al. 1998 strictly pertain only to coUisionless systems, thus 
their direct applicability to a gas-rich system like IC 2233 is somewhat unclear. 
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7.2. Position- Velocity Plots 

In Figure [11] we show a position- velocity (P-V) diagram along the major axis of IC 2233, 
together with additional P-V profiles extracted ±15" from the midplane. The major-axis 
P-V diagram of IC 2233 displays noticeable asymmetries, with the redshifted side exhibiting 
a more obvious flattening or turnover at large radii, and the blueshifted side showing an 
apparent "upturn" near x = —140" and maximum observed gas velocities ~25 km s~^ larger 
than on the redshifted side. As is frequently seen in lopsided galaxies, the "fiat" side of 
the P-V profile corresponds to the brighter horn in the global H I profile (see Figure [6] and 
§9.1j see also Swaters et al. 1999; Noordermeer et al. 2001). We also draw attention to 
the features seen in Figure [TT] near x = ±100" that indicate deviations of ~±20 km s~^ and 
—20 km s~^, respectively, from the spread of "permitted" velocities at these locations. 

It is interesting to contrast the major axis P-V plots of IC 2233 with those we derived 
for UGC 7321 in Paper I. In spite of having only a slightly smaller peak rotational velocity 
than UGC 7321, the major axis P-V plot for IC 2233 is strikingly different: while that of 
UGC 7321 is reminiscent of a "scaled-down" giant spiral galaxy (see Figure 10 of Paper I), the 
P-V profile of IC 2233 has a nearly solid-body shape, characteristic of dwarf and Magellanic 
spiral galaxies (cf. e.g., Swaters et al. 2002). At any given position along the major axis of 
IC 2233, the spread of observed velocities is narrower than what we observed in UGC 7321, 
and the peak intensities lie at smaller values of Vj-of These differences imply less rotational 
shear in IC 2233, a factor that may be important in allowing it to form the giant Hll 
complexes that are absent in UGC 7321 (see §3.3p . 

A sample of P-V diagrams extracted parallel to the minor axis of IC 2233 is shown in 
Figure [121 This figure highlights the previously mentioned high-latitude emission, extending 
to 2; ~ ±60", and the tendency for the observed velocities of this emission to lie at smaller 
values compared with the emission along the plane (see also §6.11) . This type of velocity 
structure is consistent with the expected signatures of a rotationally lagging Hi halo (e.g., 
Matthews & Wood 2003). Kinematically anomalous extraplanar gas is observed on both sides 
of the plane of IC 2233, although it is slightly more prevalent in the northeast quadrant of 
the galaxy. 

7.3. Disk Rotation Curve 

Using the data shown in Figure [TT] we have derived a major axis rotation curve for 
IC 2233 following the method described in Paper I. The result is plotted in Figure [12] and 
provided in electronic form in Table 4. The optical (Ha) rotation curve from Goad & Roberts 
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(1981) is overplotted for comparison. The generally excellent agreement between the Ha 
rotation curve and our new Hi rotation curve suggests that the latter is not significantly 
affected by beam smearing. The one notable discrepancy occurs over the interval x ~ 50"- 
100", where the Hi rotation curve is much smoother than the Ha curve. This may result 
from an irregular or patchy Hll region distribution and/or non-circular motions. Indeed, we 
see evidence of a larger spread in Hi velocities at the corresponding location (Figure [TTj) . 

A maximum rotational velocity of ~85 km s~^ is reached on both sides of the disk of 
IC 2233. Both sides also show a small but statistically significant drop in rotational velocity 
of ~10 km s~^ near the edge of the stellar disk [x ~ ±130"), just before the terminal velocity 
is reached. The locations of these dips correspond closely to where the edges of the H I disk 
are seen to warp (Figure [HD. Compared with UGC 7321, the Hi rotation curve of IC 2233 
shows a more leisurely rise relative to the stellar disk and has a less extended fiat portion. 

Despite the asymmetries visible in the major axis P-V plot of IC 2233 (Figure [TT!) and 
the global Hi profile (Figure [6]), our derived rotation curve is relatively symmetric, both 
in terms of shape and in terms of peak rotational velocity on the two sides of the disk. 
The difference in extent of the rotation curve on the two sides is also rather modest, with 
the receding side extending roughly 20" further than the approaching side. Only between 
X = 110" — 135" is a significant difference in amplitude between the two sides seen, with an 
offset of ~25 km s~^. 

A closer examination reveals that the reason why the underlying disk rotation curve 
appears fairly symmetric in spite of the asymmetries visible in Figure [11] is that the bulk 
of the underlying gas in IC 2233 has a fairly smooth and symmetric distribution (note 
the isophotes corresponding to intermediate intensities), with additional irregularities and 
asymmetries present in the form of the emission traced by the brightest, innermost contours, 
and the faintest, outermost contours. The innermost contours correspond to the bright 
clumps along the midplane visible in Figure O based on their locations in the P-V plane, 
it appears that these clumps are located at small to intermediate galactocentric radii and 
may be subject to modest non-circular motions. On the approaching side of the disk, the 
outermost contours seen in our data reveal gas extended to velocities ~25 km s~^ beyond 
the asymptotic rotational velocities shown in Figure [T31 however, this highest-velocity gas 
is of insufficient quantity to receive substantial weight in the determination of the rotation 
curve (see Paper I). 
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8. The Blue Compact Dwarf Galaxy NGC 2537 

As described in ^ the blue compact dwarf (BCD) galaxy NGC 2537 lies at a projected 
distance of 16'7 (49 dio kpc = 34 dj kpc) from IC 2233. Here we describe our new Hi 
measurements of this galaxy and briefly discuss how our observations constrain scenarios for 
its evolution (see also §9.2p . 

8.1. Optical Morphology 

NGC 2537 has been studied extensively at a variety of wavelengths, and some of its 
key properties are summarized in Table 5. This galaxy was observed at optical and near- 
infrared wavelengths by Gil de Paz et al. (2000a,2000b), who found that its nuclear region 
hosts an intermediate-age starburst {t ~30 Myr) containing numerous star-forming knots. 
Surrounding the central starburst is a diffuse, red outer disk with an underlying stellar 
population of age 5-13 Gyr, implying that NGC 2537 has been forming stars for a signiflcant 
fraction of a Hubble time. The metallicity of NGC 2537 is somewhat uncertain, but it is 
weU below solar {Z ~ 0.13 - 0.41^©; Meier et al. 2001; Wu et al. 2006). NGC 2537 has 
been imaged previously in Hi by other authors (Stil & Israel 2002a,b; Swaters et al. 2002; 
Wilcots & Prescott 2004), but the improved sensitivity and resolution of our data provide 
additional insight into the H I properties and kinematics of this galaxy. 

8.2. Properties of the Neutral Hydrogen in NGC 2537 

8.2.1. Channel Images 

In Figure [H] we present the H I channel images for NGC 2537. The same data are 
represented both as contours and as greyscale in order to highlight the existence of numerous 
small clumps of H I superposed on a lower surface density background. Many of these clumps 
are unresolved by our beam (i.e., they have sizes <500 d-j pc). The background emission 
underlying the clumps is itself rather diffuse and fragmented, and the emission within the 
individual channel images does not appear contiguous at the sensitivity limit of our data. 
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Table 5. Optical and Infrared Properties of NGC 2537 



Parameter 


Value 


Ref. 


a (J2000.0) 


08 13 14.7 


1 


5 (J2000.0) 


+45 59 26.3 


1 


Hubble type 


SB(s)m pec 


1 


Distance*^ 


5.5 - 9 Mpc 


2 


Ab (mag) 


0.232 


3 


Ar (mag) 


0.144 


3 


Measured Quantities 




1^7x1^5 


1 


ITLb^ 


12.07±0.03 


4 


B-R^ 


0.58±0.11 


4 


hr,B 


14'!4 ± 0'^4 


4 


mi(o)^ 


20.57 ± 0.09 mag arcsec"^ 


4 


Derived Quantities 




3.5 d-j kpc 


1,4 


Mb'' 


-17.16 


2,4 


T b 

^b 


8.4x10^ Lq 


2,4 




2.8x10^ dPj Lq 


1,2 



*^We adopt a distance of 7 Mpc (see §2]). 
Distance-dependent quantities are scaled in terms 
of 6^7, the actual distance in units of 7 Mpc. 

'^Corrected for foreground extinction only. 

'^Based on IRAS 60/im and lOO/^m fluxes from 
the NED database. 

References. — (1) NED database; (2) this work; 
(3) Schlegel et al. 1998; (4) Gil de Paz et al. 
2000a,b 
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8. 2. 2. The Global HI Profile of NGC 2537 

We derived a global Hi profile for NGC 2537 in the same manner as for IC 2233 ( §6.2p . 
The result is shown in Figure [T5l The uncertainty in our global profile for NGC 2537 is larger 
than that of IC 2233, both because NGC 2537 was farther from the center of the primary 
beam, and because the LSB, noncontiguous nature of its Hi emission made it more difficult 
to sum the total flux in each channel accurately. Nonetheless, consistent with the ordered 
velocity field seen in Figure [151 the derived H I profile confirms that NGC 2537 is clearly 
rotationally dominated (the measured profile width at 20% peak maximum is 121 km s~^, 
compared with the peak Hi velocity dispersion of ctv-.h: ~25 km s~^; see §8.2.5p . The Hi 
profile is somewhat lopsided, and the slopes of the two edges of the profile differ, with the 
blueshifted side having a shallower slope. The global Hi profile we have derived is in good 
agreement with the single-dish profiles published by Thuan & Martin (1981) and Huchtmeier 
& Richter (1986) and the profiles derived by Swaters et al. (2002) and Stil & Israel (2002a) 
using Westerbork aperture synthesis measurements. 

We derive an integrated Hi flux for NGC 2537 of 19.2±0.6 Jy km s^^. The uncertainty 
is a combination of a statistical error of 0.35 Jy km s^^ and an overall calibration uncertainty 
of 1%. Our integrated flux lies within the range of previously published single-dish values, 
including: 21.6±1.1 Jy km s"^ (Thuan & Martin 1981); 17.2±2.2 Jy km s^^ (Davis & 
Seaquist 1983); 20.1±1.1 Jy km s~^ (Huchtmeier & Richter 1986); 21.6±2.5 Jy km s"^ 
(Fisher & TuUy 1981); 18.5±1.2 Jy km s"^ (Tifft 1990); and 20.2±1.0 Jy km s^^ (Bottinelh 
et al. 1990). Previous aperture synthesis measurements include 18.8±0.9 Jy km s^^ (Stil & 
Israel 2002a) and 14.8±2 Jy km s'^ (Swaters et al. 2002). 

8.2.3. The Total Hi Intensity Distribution 

An H I total intensity (moment 0) image of NGC 2537 is shown in Figure dSl In the 
upper left panel, the Hi intensity contours are overplotted on a greyscale representation of 
the data, while in the upper right panel, we show the Hi contours overplotted on a blue 
image of the galaxy from the Digitized Sky Survey (DSS). Because of the very diffuse and 
patchy nature of the Hi emission in this galaxy, the percolation method described in §6.21 
picks up a large amount of low-level noise. Therefore to compute the zeroth moment we have 
summed all positive emission over the velocity range V^=368.4-502.9 km s~^. To increase 
signal-to-noise for the higher order (first and second) moments (discussed below), we rejected 
data that did not exceed 1.5 times the rms noise per channel in a version of the data that 
had been smoothed spatially with a Gaussian kernel of width 5 pixels and in velocity with 
a boxcar function of width 3 channels. 
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As seen in the upper panels of Figure [T6| we detect Hi emission in NGC 2537 extending 
to ~1.5-2 times the hmiting isophote on the blue DSS image, with the extent of the Hi 
relative to the starlight varying with azimuth. The mean ratio of the H I diameter to the D25 
optical diameter is ~3.5. As is evident from the individual channel images (Figure [T^ . the 
global morphology of the neutral gas is quite clumpy and fragmented, and several patches 
within the disk appear nearly devoid of Hi emission. This flocculent appearance of the 
Hi is not confined to the star-forming portion of the galaxy, implying that it is not linked 
to feedback processes or to the instabilities giving rise to star formation. The overall Hi 
distribution and morphology of NGC 2537 are similar to other low-mass starburst and post- 
starburst galaxies observed with comparable spatial resolution (e.g., Israel & van Driel 1990; 
Meurer et al. 1996; van Zee et al. 1998; Wilcots & Miller 1998; Hunter et al. 1999; Thuan 
et al. 2004), although frequently BCDs show global Hi morphologies that are much more 
disorganized and irregular than seen in NGC 2537 (cf. van Zee et al. 1998,2001; Thuan et 
al. 2004). 

A notable feature of the outer H I disk of NGC 2537 is the spiral arm-like feature on 
the eastern edge, which extends nearly 180° around the galaxy. This feature has no optical 
counterpart, and as seen in §8.2.41 the gas in this arm follows the same ordered rotation 
as the rest of the disk. One possibility is that this feature could be a remnant from a past 
interaction ( §9.ip . However, gaseous arms induced by interactions tend to be short-lived, 
persisting only one or two rotation periods (e.g., Quinn 1987; lono et al. 2004), implying 
that the interloper should still be nearby. As discussed in §2] and §9.11 we are unable to 
identify a likely culprit. Moreover, the tight winding of the NGC 2537 arm, together with 
its regular rotation and lack of a counter-arm all are inconsistent with classic tidal features 
(see Toomre & Toomre 1972). Explaining the arm through more indirect dynamical effects, 
such as swing amplification (Toomre 1981), also appears problematic given the extremely 
low H I surface density of the disk. The Toomre Q parameter of the outer disk is quite high 
{Q >10) implying a high degree of stability against axisymmetric perturbations; in such 
a case, the swing amplification mechanism will be very inefficient, although this problem 
could be overcome if the H I disk contains sufficient dark matter to render it gravitationally 
unstable (see the discussion in Masset & Bureau 2003). 

The most striking feature of the inner Hi disk of NGC 2537 is the bright "ring" 
approximately 1' (2.0^7 kpc) in diameter. Typical column densities along the ring are 
nni ~ 7.6 x 10^*^ cm~^, although several brighter knots with nni ~ 9.5 x 10^*^ cm~^ are 
superposed. In contrast, inside the ring the mean column density is roughly a factor of two 
lower (~ 4.7 x 10^*^ cm~^), although part of the gas inside is likely to be molecular. One of 
the regions of bright optical continuum emission in NGC 2537 lies coincident with the center 
of the Hi ring, but based on the Ha and [Olll] images of Gil de Paz et al. (2000b), there 
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does not appear to be any significant amount of ionized gas at tliis location. 

Using tfie data from Figure [TBk . we liave derived a radial H I surface density profile 
for NGC 2537 by taking measurements within a series of concentric, elliptical annuli. The 
center of the annuli was taken to be the kinematic center listed in Table 6; this position 
approximately corresponds to the center of the bright H I ring in Figure [T6h.b. The ellipticity 
of the annuli was taken to be 0.1 and their major axis position angle to be 172°.0 (see §8.2.41 
below). The results are presented in Figure [T7i We have not corrected the inferred Hi 
surface density for inclination since this correction is small (~18%). We see that the Hi 
surface density in the outer disk of NGC 2537 is extremely low, comparable to values found 
in typical dwarf irregular and LSB spiral galaxies, including IC 2233 (Figure [8]). However, in 
contrast to IC 2233, NGC 2537 shows a steep rise in Hi density toward its center, where the 
Hi surface density reaches a peak of ~ IIMq pc~^ — comparable to values seen in the inner 
parts of normal late- type spiral galaxies (Cayatte et al. 1994). Similar enhancements in the 
central H I surface density are also observed in other BCD systems (Meurer et al. 1996; van 
Zee et al. 1998). 

We have overplotted as a dashed line on Figure [T71 the canonical star formation threshold 
of Kennicutt (1989). We have assumed that the gas is isothermal with a constant sound 
speed of 7 km s~^. We see that the computed Kennicutt threshold lies significantly above 
the measured gas surface density at all radii. However, scaling the measured Hi curve 
by a factor of 1.34 to account for He and by an additional factor of 1.9 reconciles it with 
the computed instability threshold over the portion of the galaxy where stars are visible 
(delineated by the vertical bar). The latter scaling factor is comparable to the expected 
contribution to the total gas surface density from the molecular gas: Gil de Paz et al. (2002) 
measured ~ (0.4-5) x 10^ Mq within the the central ~ 0.7^7 kpc (r <20") of NGC 2537. 
Therefore, in contrast to IC 2233, the locations of star formation within NGC 2537 appear 
to be in reasonable agreement with Kennicutt 's star formation criterion so long as the total 
gas surface density is considered. 



8.2.4. The Hi Kinematics of NGC 2537 

The H I velocity field shown in Figure [16] reveals that NGC 2537 exhibits well-ordered 
rotation. Note that the H I arm mentioned in the previous section participates in the regular 
disk rotation (although a series of kinks can be seen in the isovelocity contours as they 
cross the region between the arm and the main disk). While BCDs as a class are in general 
rotationally-supported, only a small fraction exhibit the type of well-ordered velocity field 
seen in NGC 2537 (cf. van Zee et al. 1998,2001; Thuan et al. 2004). 
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Table 6. 21-cm Radio Properties of NGC 2537 



Parameter 



Value 



Measured Quantities 



a (J2000.0)'' 
(5 {J2000.0)'' 

Peak Hi column density'' 
Dm 

W20 
Wso 

Fcont (21-cm) 
Position angle° 
Inclination'' 

V ^ 

V max 



08 13 14.84 (±0.04) 
+45 59 30.5 (±0.4) 
9.5 X 10^° atoms cm~^ 
5.'90±0.'05 
19.2 ±0.6 Jy km s"' 
120.9 km s-i 
99.2 km s~i 
445.3 ±0.2 km s-^ 
443.5 km s"! 
10.3 ±0.6 mjy 
172.0° ± 0.2° 
33° ± 1° 
87 km s-i 



"Mh? 

Mm /Lb 
Dm/D25 



Derived Quantities 



(2.22 ± 0.07) X 10** 4 ^0 
9.5x10^ d7 Mq 
2.5 Mq/Lq 
43 ^ 
3.5 



Trom tilted-ring rotation curve fit. 

''At 16" resolution. 

^Heliocentric; optical definition. 

''From global Hi profile (optical definition). 

"Based on the kinematic major axis of the Hi velocity 
field (see ^STZit 

'Assuming the H I is optically thin and dj is the distance 
to NGC 2537 expressed in units of 7 Mpc. 

gProm A/dyn = 2.326 X lQi^rV'^(r), where we have taken 
r = 5.4 d7 kpc (160") and V{t) = Vmax =87 km s"'. 
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Table 7. NGC 2537 Rotation Curve 



IVIajor 3:xis distanc© (arcsec) 


Circu.la.r Velocity (km s 


Uiicerta.iiity (km s ^) 


-165.0 


497.5 


3.2 


-156.0 


492.0 


2.0 


-147.0 


491.2 


3.2 


-138.0 


493.5 


3.1 


-129.0 


485.8 


5.3 


-120.0 


603.2 


7.6 


-111.0 


497.8 


5.4 


-102.0 


493.2 


2.8 


-93.0 


487.8 


3.8 


-84.0 


486.9 


1.8 


-75.0 


488.6 


1.1 


-66.0 


492.9 


0.7 


-57.0 


490.6 


0.9 


-48.0 


486.6 


0.8 


-39.0 


484.4 


0.8 


-30.0 


483.8 


0.9 


-21.0 


479.2 


1.0 


-12.0 


471.8 


1.4 


-3.0 


461.8 


2.4 


6.0 


442.4 


2.8 


15.0 


417.4 


2.1 


24.0 


410.0 


1.3 


33.0 


411.4 


2.8 


42.0 


411.4 


8.0 


51.0 


412.5 


5.2 


60.0 


409.4 


1.9 


69.0 


400.0 


3.0 


78.0 


402.6 


1.3 


87.00 


403.1 


1.9 


96.0 


394.6 


5.6 


105.0 


384.6 


8.0 


114.0 


408.0 


1.3 


123.0 


396.0 


5.3 


132.0 


401.1 


1.5 


141.0 


404.0 


1.7 


150.0 


405.1 


1.5 


159.0 


395.1 


1.6 
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In Figure [18] we show a P-V diagram extracted along the kinematic major axis of 
NGC 2537 (lying at PAf» 172°; see below). To improve signal-to- noise, we have averaged over 
a 15"-wide strip. Using the data shown in Figure [T5| we have also measured a disk rotation 
curve for NGC 2537 by fitting intensity-weighted Gaussians to a series of one-dimensional 
slices extracted at 9" intervals along the major axis. The result is shown in Figure [19] and 
provided in electronic form in Table 7. Despite the rather high central velocity dispersion 
observed in NGC 2537 (see below) we have not applied corrections to the observed rotation 
curve for asymmetric drift. We estimate that such corrections are likely to be significant 
(~ 20 — 30%) only within the centralmost regions of the galaxy (r <30"), and the corrections 
to individual points are difficult to compute accurately with only a few resolution elements 
across this portion of the galaxy (cf. Meurer et al. 1996; Masset & Bureau 2003). Even 
without correcting for asymmetric drift, we find the rotation curve of NGC 2537 rises rather 
steeply, reaching a plateau within ~ 30" (~ 1.0 kpc, i.e., two beam widths). Beyond this, 
a more extended region is seen on both sides of the disk. This extended portion is nearly 
flat, although it continues to rise slightly with increasing galactocentric distance. Superposed 
on the outer, nearly flat portion of the rotation curve, a series of "wiggles" with amplitude 
~10 km s^^ is visible. The onset of these velocity undulations corresponds with the observed 
edge of the stellar disk. 

Based on the rotation curve in Figure [T9l the peak obs ervediotaiion velocity of NGC 2537 
is ~50 km s^^. Correcting this to the true rotational velocity requires an estimate of the 
inclination of the galaxy. The morphology of this BCD makes an inclination determination 
from optical images difficult, as the intrinsic flattening of this type of galaxy is poorly known. 
To better constrain the inclination and other kinematic properties of NGC 2537, we have 
used the velocity-field analysis method of van Moorsel & Wells (1985) which models the 
galaxy as a rotating thin disk. Assuming an exponential parameterization of the rotation 
curve leads to rapid convergence but also shows that the bright central ring is tilted with 
respect to the rest of the galaxy. Fitting the region outside the central 40" (in radius) re- 
sults in an inclination angle of 33 ± 1°, whereas the region of bright Hi clumps is well fit 
by a rotating ring with an inclination of 28 ± 3°. Both regions yield a consistent center of 
rotation at ajsooo.o = 08^* 13"^ 14.84^ (±0.04^); 5j2ooo.o = +45° 59' 30"5 (±0':4), where the 
errors are formal values derived from the fits. The systemic velocities derived from both fits 
are consistent as well, yielding an average value of 445.3 ± 0.2 km s~^. The major axis is 
found to lie at PA= 172.0° ± 0.2°, with both fits again showing good agreement. Finally, 
the deprojected rotational velocities are 80 ± 7 km s~^ for the "ring" and 87 ± 2 km s~^ for 
the asymptotic (maximum) rotational velocity of the disk of the galaxy. The two velocity 
curves join smoothly at the radial distance of 40". This value of V^max is surprisingly high 
for an optically tiny galaxy like NGC 2537, and it exceeds most of the values measured for 
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other BCDs, which tend to be a factor of two smaller (van Zee et al. 1998), although a few 
exceptions are known (e.g., Meurer et al. 1996). Nonetheless, the rather large derived 
for NGC 2537 is congruous with the shape of its rotation curve, which shows a turnover and 
extended flat portion as opposed to the solid body-type rotation curve seen in lower-mass 
dwarfs. 



8.2.5. HI Velocity Dispersion 

A map of the line-of-sight H I velocity dispersion in NGC 2537 is shown in Figure [T6l i. 
A radial gradient in the velocity dispersion is apparent, with values ranging from ay hi ~4- 
8 km s~^ in the galaxy outskirts to as high as crv,m ~ 25 km s~^ near the center. The values 
in the outskirts are comparable to the canonical values of ~6-9 km s~^ typical of both dwarf 
and spiral galaxies; however, the values near the disk center reach a signiflcant fraction of 
the rotational velocity in the disk, implying appreciable pressure support. The gas with the 
highest measured velocity dispersions coincides with the location of the bright star forming 
regions in the stellar disk of NGC 2537, although a few areas with av,m > 10 km s~^ are 
also seen outside the stellar disk. A similar pattern of H I velocity dispersion has also been 
observed in other BCDs, including NGC 2915 (Meurer et al. 1996). In this latter case, 
Wada et al. (2002) argued that a stellar bar is responsible for driving the high central 
turbulence. However, NGC 2537 does not show any sign of a stellar bar, implying that 
another mechanism (most likely feedback from star formation) is responsible for the elevated 
velocity dispersion. 

9. Discussion 

9.1. The Origin of the Structural and Kinematic Features of the IC 2233 Disk 

In earlier sections we have drawn attention to a variety of interesting features of the 
Hi distribution and kinematics of IC 2233, including warping and flaring of the Hi layer; a 
corrugated vertical structure; lopsidedness in the H I intensity distribution, global H I proflle 
and P-V diagrams; and evidence for vertically-extended, rotationally anomalous gas. Here 
we discuss possible origins for these various features and note the likelihood that some or all 
of these phenomena may share a common origin. 

Some of the features of IC 2233 are known to be quite common in disk galaxies. For 
example, the study of Garcia- Ruiz et al. (2002) established that warping is nearly ubiquitous 
in disk galaxies and affects nearly 100% of galaxies where the Hi disk extends beyond 
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the stellar distribution. Similarly, lopsidedness in the stellar and gas distributions is now 
well known to affect a significant fraction of galaxies (Richter & Sancisi 1994; Zaritsky & 
Rix 1997; Haynes et al. 1998), and asymmetries seem to be even more common among 
extreme late-type disks, affecting more than half of such systems (Matthews et al. 1998). 
To date, rotationally anomalous extraplanar Hi and disk corrugations have been confirmed 
in only a handful of galaxies (see §6.4[§7.2p . but relatively few investigators have searched 
systematically for these phenomena, thus their frequency of occurrence is not yet known. 

Tidal interactions and minor mergers frequently have been invoked to explain a variety 
of structural and kinematic phenomena in galaxies. For example, tidal interactions have been 
suspected as a key trigger of galaxy lopsidedness (e.g., Zaritsky & Rix 1997) and have been 
proposed as a possible driver of warps (e.g., Weinberg 1995) and of disk corrugations (e.g., 
Edelsohn & Elmegreen 1997). This raises the possibility that this process could be shaping 
IC 2233. Indeed, we chose to observe IC 2233 as a suspected "victim" of a tidal encounter 
(see [H]). However, numerical simulations predict that most tidally-induced disturbances 
will be rather short-lived, and prevalent only during the interval when the minor merger is 
ongoing (Walker et al. 1996; Sellwood et al. 1998; Noguchi 2001; Bournaud et al. 2005). 
Our observations do not support this scenario, as we find no neighbor with sufficient mass 
to induce tidal effects in IC 2233, and we have not identified any stellar or gaseous remnant 
of a possible intruder. 

In ^ we noted the lack of any candidates for bona fide companions to IC 2233 based 
on optical data and argued that its previously suspected neighbor, NGC 2537, most likely 
lies at a different distance. Furthermore, our Hi observations of both galaxies show that 
even if they lay at similar distances, NGC 2537 and IC 2233 would not constitute a bound 
pair. Assuming a minimum separation (equal to the projected separation at 10 Mpc), a 
relative velocity equal to the difference in the recessional velocities, and total galaxy masses 
equal to the dynamical masses given in Tables 3 and 6, then the kinetic energy of the system 
(~ 7.2 X 10^^ erg) would exceed the potential energy (~ 2.6 x 10^^ erg) by roughly a factor of 3. 
Furthermore, the velocity separation of the centers of mass of the galaxies exceeds their "flat" 
rotational velocities, again indicating that the galaxies are not gravitationally bound. Finally, 
a comparison of Figures [TD] and [12] shows that any interaction would have been retrograde, 
thus minimizing its dynamical effects (e.g., Velazquez & White 1999). Consistent with these 
arguments, our deep Hi observations of IC 2233 and NGC 2537 have failed to uncover any 
overt signatures of a recent or ongoing interaction, such as tidal debris or counter-rotating 
gas. 

To investigate the environment of IC 2233 further, we have also performed a systematic 
search for gas-rich neighbors using the matched-filter method described in Paper I. We 
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explored a 30' x 30' region over the velocity interval covered by our observations (see Table 2) 
and searched for signals with velocity widths of up to ~40 km s~^. This search uncovered 
seven features with amplitudes > 4a, but none greater than 5a. Four of the Aa features 
are positive (emission) features, while three are negative, suggesting these simply represent 
the expected Gaussian tail of the noise distribution. Moreover, none of the positive features 
have optical counterparts in the SDSS data (see We therefore do not consider any of 
these to be likely candidates for companions to IC 2233. Based on these results, we derive 
an {8a) upper limit to the Hi mass of any companions to IC 2233 of Mhi< 3 x lO^df^ Mq 
within a radius of 30' (~ 90 rfio kpc). Similar limits for NGC 2537 are Mhi< 1.5 x lO^d^ Mq 
within ~ 60 dj kpc. 

If the properties of IC 2233 cannot readily be accounted for by a recent interaction, it 
is necessary to consider alternative explanations. For example, Levine & Sparke (1998) and 
Noordermeer et al. (2001) have proposed that galaxy lopsidedness may be linked to the disk 
lying offset from the center of its global potential. Noordermeer et al. (2001) presented some 
predictions of this model that can be compared with observations. The details depend on 
whether the disk rotates prograde or retrograde relative to the halo. However, in general the 
models predict: (1) modest asymmetries in shape and extent of the two sides of the major 
axis P-V curve, with the more extended side showing a slight velocity decrease and the less 
extended side showing a slight "upturn" near the last measured point; (2) a global Hi profile 
with a higher peak on the side corresponding to the more extended side of the rotation curve 
and a lower, sloping horn on the opposite side; (3) various degrees of twisting or distortion 
of the velocity field (which depend on the orientation of the axis of symmetry of the disk 
relative to the observer). All of these signatures are observed in IC 2233 (Figures [6], [T0| 
[TTIl . Unfortunately, Noordermeer et al. (2001) did not offer specific suggestions as to what 
mechanism might lead to an offset between the disk and the overall galactic potential. 

There are intriguing similarities between the velocity field, rotation curve, and global 
Hi profile of IC 2233 and those of the dwarf LSB spiral NGC 4395 (see Swaters et al. 
1999; Noordermeer et al. 2001), a galaxy we have suggested might be an analog of IC 2233 
observed face-on (§6.3). NGC 4395 is not known to be tidally interacting or undergoing 
a minor merger, hence these similarities suggest the possibility of a more universal and 
continuous driver for these traits in low-mass, late-type disks (see also Wilcots & Prescott 
2004). Spontaneously arising asymmetries (e.g., those linked with spiral arm formation) 
are one possibility, although Bournaud et al. (2005) argue that in most instances these are 
likely to be too weak; these authors suggest instead that a more promising explanation is 
the continuing cosmological accretion of gas. Cosmological models predict that such infall 
should continue to the present day, particularly in low-density field environments such as 
those inhabited by IC 2233 and NGC 4395 (e.g., Keres et al. 2005). 
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The possibility that IC 2233 might be continuing to accrete material from the IGM is 
intriguing in light of several other recent theoretical studies that predict that cosmological 
accretion can account for a variety of properties observed in isolated galaxies, including the 
presence of anomalously rotating extraplanar gas (Fraternali et al. 2007) and warping of the 
disk (Jiang & Binney 1999; Shen & Sellwood 2006). While none of these studies specifically 
investigated whether disk corrugations might also be triggered by slow gas accretion, pre- 
vious studies linking corrugations with warps (Nelson 1976; Sparke 1995; Masset & Tagger 
1997) and velocity undulations with disk lopsidedness (Schoenmakers et al. 1997) suggest 
this would be an interesting topic for further research. Additional high-resolution, high- 
sensitivity Hi studies of apparently isolated, low-mass spirals would also help to establish 
how ubiquitous the simultaneous appearance of these various phenomena is in the absence of 
close neighbors and permit more rigorous comparisons with the morphological and kinematic 
signatures of gaseous infall predicted by models. 

9.2. A Comparison between NGC 2537 and IC 2233: Clues to their 

Evolutionary Histories 

As emphasized by Matthews & Gallagher (1997), one intriguing aspect of extreme late- 
type disk galaxies is the existence of such a diversity of disk morphologies within a relatively 
narrow range of mass and luminosity (see also Noguchi 2001). IC 2233 and NGC 2537 are 
two excellent illustrations of this contrast; both are rotationally-dominated galaxies having 
similar peak rotational velocities, dynamical masses, and blue luminosities (see Tables 3 & 
6), yet they exhibit a variety of differences in terms of their optical and Hi morphologies. 
We now briefly comment on the possible origin and implications of these differences. 

One key difference between IC 2233 and NGC 2537 becomes apparent from a compar- 
ison of their respective disk rotation curves (Figures [T3] fc [T9l) . Compared with IC 2233, 
the rotation curve of NGC 2537 rises far more steeply in the central regions, implying a 
much higher central mass concentration. Indeed, such enhanced central mass concentrations 
appear to be a hallmark of BCDs (e.g., van Zee et al. 1998). 

Previous authors have suggested that BCDs might represent a brief starbursting state 
of ordinary dwarf galaxies (e.g., Taylor et al. 1993; Papaderos et al. 1996). While we cannot 
rule out such a scenario for some low-mass BCDs, our present observations suggest that it 
is unlikely that a galaxy like NGC 2537 will ever evolve into a "normal" LSB disk galaxy 
(e.g., similar to IC 2233), or vice versa. 

If a tidal interaction or minor merger were invoked to explain a transition from LSB disk 
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to starbursting BCD (e.g., Taylor 1997), one difficulty would be to account for the properties 
of the outer gas disk of NGC 2537, which appears quite dynamically cold with no signs of 
recent dynamical heating or perturbations. Moreover, in contrast to typical gas-rich LSB 
spirals — which tend to have blue outer disks, shallow light profiles, and little central light 
concentration — the underlying stellar disk of NGC 2537 is old, red, and compact; thus it is 
unclear that the required pre-merger progenitor would resemble any ordinary class of disk 
galaxy. Finally, we note that NGC 2537 seems to join an increasing number of BCDs that 
are found to reside in rather isolated environments (Telles & Terlevich 1995; van Zee et al. 
1998). 

As an alternative means of linking BCD and LSB/dwarf disk galaxies, Papaderos et 
al. (1996) proposed that the reverse transition (from BCD to LSB) might occur via secular 
evolution processes. However, this would necessitate a large-scale redistribution of mass, 
presumably via starburst-driven winds followed by cooling and infall of this gas from the halo. 
While this might be possible for extremely low-mass dwarfs, this requirement is particularly 
problematic for NGC 2537, since its gravitational potential is too steep to permit substantial 
mass-loss via mechanical energy from stellar winds and supernovae (see e.g., Ferrara & 
Tolstoy 2000). 

Given the above difficulties, the current body of observations of NGC 2537 seems to 
be most consistent with the suggestion that at least some BCDs are a special subset of 
low-mass disk galaxies whose central mass concentrations lie at the extreme end of the 
present-day distribution (Salzer & Norton 1999; Ferrara & Tolstoy 2000; Noguchi 2001). 
Models by Noguchi (2001) predict that viscous evolution would be particularly efficient in 
such centrally-concentrated galaxies, thus giving a natural explanation for their ability to 
funnel gas inwards to fuel episodic central starbursts; in contrast, his models predict that 
the low central mass densities of LSB galaxies like IC 2233 will help them to maintain a 
fairly flat and largely subcritical H I surface density across much of their disk. We conclude 
that IC 2233 and NGC 2537 are therefore examples of small disk galaxies that have evolved 
in relative isolation, and whose inherent structural differences have allowed them to preserve 
important differences in their morphologies and star-forming properties. 

10. Summciry 

We have presented VLA H I imaging of the cdgc-on Sd 'superthin' spiral galaxy IC 2233 
as well as the blue compact dwarf (BCD) NGC 2537. We have also described results from 
new optical B, R, and Ho; imaging and photometry for IC 2233. 
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We confirm that IC 2233 is an intrinsically LSB galaxy, having a deprojected cen- 
tral surface brightness yUB,j(0) ~ 22.6 mag arcsec"^ and an extremely low Hi surface den- 
sity over its entire disk (Shi <3M0 pc~^). Similar to many other late-type, LSB galax- 
ies, the Hi component of IC 2233 comprises a significant fraction of its observed baryons 
(Mhi/Lb ~0.83 Mq/Lq). Despite evidence for localized ongoing star formation in IC 2233 
in the form of prominent Hll complexes and shells, red and blue supergiant stars, and a 
very blue integrated color [{B — R)q = 0.67 ± 0.15], we detect no significant 21-cm radio 
continuum emission from the galaxy, and the global star formation rate inferred from the Ha 
emission and from the IRAS FIR emission both imply an extremely low globally-averaged 
star formation rate ( <0.05 Mq yr~^). 

Both the Hi and ionized gas disks of IC 2233 are clumpy and vertically distended, 
exhibiting scale heights comparable to that of the young stellar disk. The thickness of both 
the gas and the stars flares with increasing galactocentric radius. We have quantified the 
flaring of the Hi disk and estimate an increase in scale height of a factor of two across the 
galaxy. We find evidence that IC 2233 also contains a component of "anomalous" extraplanar 
Hi emission whose rotation does not follow that of the material in the midplane. Future 
three-dimensional kinematic modeling should help to confirm whether this material comprises 
a part of a rotationally lagging H I halo similar to those now confirmed in a number of other 
spiral galaxies. 

The H I disk of IC 2233 exhibits a mild lopsidedness as evidenced by differences in the 
Hi intensity distribution on the two sides of the disk, differences in shape and extent of 
the major axis P-V curves on the two sides of the galaxy, and by the asymmetric shape 
of the global Hi profile. The origin of this lopsidedness is unclear, although its kinematic 
signatures are consistent with model predictions for a disk lying offset from the center of its 
overall halo potential. 

A particularly intriguing feature of the H l disk of IC 2233 is the presence of a global 
corrugation pattern with a period of ~ 7 dio kpc and an amplitude of ~ 150 dio pc. These 
undulations may represent bending instabilities or be linked to underlying spiral structure. 
Their presence suggests that the optically diffuse disk of IC 2233 is largely self-gravitating. 
Outside of the stellar disk, a mild warp of the H I disk is also observed. 

The properties of IC 2233 form an interesting contrast to those of the LSB superthin 
galaxy UGC 7321 (Paper I). Although IC 2233 is only a factor of two less massive than 
UGC 7321, its kinematic characteristics resemble those of a solid-body rotator ( §7.2p . while 
those of UGC 7321 are more similar to those of a giant spiral. The stellar disk as well as the 
neutral and ionized gas layers of IC 2233 are intrinsically thicker than those of UGC 7321 
and exhibit far more complex structure. These differences may reflect in part the weaker 
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self-gravity of the IC 2233 disk and the decreasing importance of rotation shear. The re- 
spective properties of IC 2233 and UGC 7321 are consistent with previous suggestions that 
important changes occur in the structural and ISM properties of galaxies over the mass 
interval corresponding to V^ot ~100 km s~^ (Dalcanton et al. 2004; Matthews et al. 2005). 

We have confirmed that the BCD NGC 2537 is a rotationally-dominated galaxy and 
estimate its maximum (dcprojcctcd) rotational velocity to be 14nax ~87 km s"-*". This is at 
the high end of values observed for BCDs. The steep inner rise of its rotation curve indicates 
a high central matter density. 

The H I velocity field of NGC 2537 is quite regular and shows no obvious signs of recent 
perturbations. An arm-like feature in the outer H I disk shows regular rotation and does not 
appear to be tidal in origin. The Hi disk of NGC 2537 extends to ~3.5 times the (-D25) 
diameter of its stellar disk, and its outer regions appear dynamically cold and fragmented. 
In the inner Hi disk we measure Hi velocity dispersions as high as av,m ~ 25 km s~^, 
indicative of significant turbulent motions. The inner disk also exhibits a tilted H I ring with 
a diameter of ~ 2.0 dj kpc and peak column density of ~ 10^^ atoms cm~^ — several times 
higher than that of the surrounding gas. 

Several of the properties of of IC 2233, including its lopsidedness, warped H I disk, disk 
corrugations, and the presence of anomalously rotating, vertically extended Hi emission, 
might seem readily explained by a recent interaction or minor merger, but we have not 
uncovered any direct evidence of such an event. We also argue that NGC 2537 is unhkely to 
be the required perpetrator. Indeed, both IC 2233 and NGC 2537 each appear to be evolving 
in relative isolation. Although IC 2233 and NGC 2537 have sometimes been described as 
a "binary pair" in the past literature, recent distance estimates suggest they probably lie 
at different distances. Moreover, even if they lay at the same distance, our data show they 
would not form a bound pair and any encounter would have been retrograde, implying that 
any mutual interaction would have had minimal impact on the properties of the two galaxies. 
Our search for other companions to both galaxies using imaging and spectroscopic data from 
the SDSS as well as our new H I data have failed to identify any additional candidates for close 
neighbors. Slow, continuing accretion of intergalactic gas may be one means of accounting 
for a number of the observed properties of IC 2233 in the absence of a companion. In the 
case of NGC 2537, the intrinsically high central matter density of the galaxy coupled with 
efficient viscous evolution may account for its bursting state without the need for a recent 
interaction. 

We thank J. Gallagher for assistance with the WIYN observations and A. Seth for 
computing the distance modulus to IC 2233 from his HST data. We also acknowledge 



-48- 



valuable discussions with a number of colleagues, including D. Hogg, J. Hibbard, M. Roberts, 
L. Sparke, T. Cornwell, J. Condon, B. Cotton, and E. Greisen. LDM was supported by a 
Jansky Fellowship from NRAO and a Clay Fellowship from the Harvard-Smithsonian CfA. 
The VLA observations described here were part of program AM649. This research has 
made use of: the NASA/IPAC Infrared Science Archive and the NASA/IPAC Extragalactic 
Database (NED) which are operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and Space Administration, 
and the Digitized Sky Surveys (DSS), which were produced at the Space Telescope Science 
Institute under U.S. Government grant NAG W-2166. Finally, we thank the referee for a 
thoughtful and thorough reading of our paper. 



REFERENCES 

Abadi, M. G., Navarro, J. F., Steinmetz, M., Eke, V. R. 2003, ApJ, 591, 499 
Alfaro, E. J. & Efremov, Y. N. 1996, RMxAC, 4, 1 

Alfaro, E. J., Perez, E., Gonzalez Delgado, R. M., Martos, M. A., & Franco, J. 2001, ApJ, 
550, 253 

Alton, P. B., Davies, J. I., & Bianchi, S. 1999, A&A, 343, 51 
Baldwin, J. E., Lynden-Bell, D., & Sancisi, R. 1980, MNRAS, 193, 313 
Bell, E. F. 2003, ApJ, 586, 794 

Bergvall, N. & Ronnback, J. 1995, MNRAS, 273, 603 

Binney, J. & Tremaine, S. 1987, Galactic Dynamics, (Princeton: Princeton University Press) 
Bizyaev, D. & Mitronova, S. 2002, A&A, 389, 795 

Bos, A. 1984, in Indirect Imaging, ed. J. A. Roberts, (Cambridge: Cambridge University 
Press), 239 

Bos, A. 1985, PhD dissertation. University of Groningen 

Bosma, A. 2002, in Disks of Galaxies: Kinematics, Dynamics, and Perturbations, ed. E. 
Athanassoula, A. Bosma, and R. Mujica, ASP Conf. Series, Vol. 275 (San Fran- 
cisco: ASP), 23 

Bottinelli, L., Gougucnhcim, L., Fouque, P. & Paturel, G. 1990, A&AS, 82, 391 

Bournaud, F., Combes, F., Jog, C. J., & Puerari, I. 2005, A&A, 438, 507 

Bridges, T. J., Carter, D., Harris, W. E., & Pritchet, C. J. 1996, MNRAS, 281, 1290 



- 49 - 



Briggs, D. S. 1995, Ph.D. Dissertation, New Mexico Institute of Mining and Technology, 
Socorro, NM (http://www.aoc.nrao.edu/ftp/dissertations/dbriggs/disslitmII), p. 64 

Brinks, E., Walter, F., & Ott, J. 2002, in Disks of Galaxies: Kinematics, Dynamics, and 
Perturbations, ASP Conf. Series, Vol. 275, ed. E. Athanassoula, A. Bosma, & R. 
Mujica, 57 

Burton, W. B. 1988, in Galactic and Extragalactic Radio Astronomy, Second Edition, ed. 
G. L. Vershuur and K. I. Kellermann (New York: Springer- Verlag) , 295 

Cayatte, V., Kotanyi, C., Balkowski, C., & van Gorkom, J. H. 1994, A J, 107, 1003 

Condon, J. J. 1992, ARA&A, 30, 575 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. P., Perley, R. A., Taylor, G. B., & 
Broderick, J. J. 1998, AJ, 115, 1693 

Cornwell, T. J., Uson, J. M., & Haddad, N. 1992, A&A, 258, 583 

Dalcanton, J. J., Yoachim, P., & Bernstein, R. A. 2004, ApJ, 608, 189 

Davis, L. E. & Seaquist, E. R. 1983, ApJS, 53, 269 

de Vaucouleurs, G., de Vaucouleurs, A., & Corwin, H. G. 1976, Second Reference Catalogue 
of Bright Galaxies (Austin: University of Texas Press) 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G. Jr., Buta, R. J., Paturel, G., & 
Fouque, P. 1991, The Third Reference Catalogue of Bright Galaxies (New York: 
Springer) 

Dickey, J. M. 1996, in Unsolved Problems of the Milky Way, lAU Symp. 169, ed. L. Blitz 
& P. Teuben, 489 

Dohm-Palmer, R. C. & Skillman, E. D. 2002, AJ, 123, 1433 

D'Onghia, E. & Burkert, A. 2004, ApJ, 612, 13 

Edelsohn, D. J. & Elmegreen, B. G. 1997 MNRAS, 287, 947 

Eliche-Moral, M. C, Balcells, M., Aguerri, J. A. L., & Gonzalez- G aria, A. C. 2006, A&A, 
457, 91 

Elmegreen, B. G. & Hunter, D. A. 2000, ApJ, 540, 814 
Elmegreen, B. G. & Parravano, A. 1994, ApJ, 435, L121 
Ferguson, A. M. N. 1998, PhD thesis, Johns Hopkins University 

Ferguson, A. M. N., Wyse, R. F. G., Gallagher, J. S., & Hunter, D. A. 1998, ApJ, 506, L19 
Ferrara, A. & Tolstoy, E. 2000, MNRAS, 313, 291 
Fisher, J. R. & Tully, R. B. 1981, ApJS, 47, 139 



-50- 



Fleurier, C. & Chapelle, J. 1974, Computer Physics Communications, 7, 200 

Florido, E., Battaner, E., Gros, A., Prieto, M., & Mediavilla, E. 1992, Ap&SS, 190, 293 

Florido, E., Battaner, E., Sanchez-Saavedra, M. L., Prieto, M., & Mediavilla, E. 1991, 
MNRAS, 251, 193 

Fraternali, F., Binney, J., Oosterloo, T., Sancisi, R. 2007, NewAR, 51, 95 

Fraternali, F., Oosterloo, T., Sancisi, R., & van Moorsel, G. A. 2001, ApJ, L47 

Fuchs, B. 2002, in Dark Matter in Astro- and Particle Physics, ed. H. V. Klapdor-Kleingrothaus 
and R. D. Viollier, (Berlin: Spriner), 28 

Fuchs, B. 2003 Ap&SS, 284, 719 

Garcia-Ruiz, I., Sancisi, R., & Kuijken, K. 2002, A&A, 394, 769 
Garrido, O. Marcehn, M. & Amram, P. 2004, MNRAS, 349, 225 
Gerritsen, J. P. E. & de Blok, W. J. G. 1999, A&A, 342, 655 

Gil de Paz, A., Silich, S. A., Madore, B. F., Sanchez Contreras, C, Zamorano, J., & Gallego, 
J. 2002, ApJ, 573, LlOl 

Gil de Paz, A., Zamorano, J., & Gallego, J. 2000a, A&A, 361, 465 

Gil de Paz, A., Zamorano, J., Gallego, J., & Dommguez, F. de B. 2000b, A&AS, 145, 377 
Goad, J. W. & Roberts, M. S. 1981, ApJ, 250, 79 
Griv, E. & Chiueh, T. 1998, ApJ, 503, 186 

Gum, C. S., Kerr, F. J., & Westerhout, G. 1960, MNRAS, 121, 132 

Haynes, M. P., Giovanelli, R., Chamaraux, P., da Costa, L. N., Freudhng, W., Salzer, J. J., 
& Wegner, G. 1999, AJ, 117, 2039 

Haynes, M. P., van Zee, L., Hogg, D. E., Roberts, M. S., & Maddalena, R. J. 1998, AJ, 115, 
62 

Hofner, P. & Sparke, L. S. 1994, ApJ, 428, 466 

Huchtmeier, W. K. & Richter, O.-G. 1986, A&AS, 63, 323 

Hunter, D. A., van Woerdcn, H., & Gallagher, J. S. 1999, AJ, 118, 2184 

lono, D., Yun, M. S., & Mihos, J. C. 2004, ApJ, 616, 199 

Israel, F. P. & van Driel, W. 1990, A&A, 236, 323 

Jiang, I. G. & Binney, J. 1999, MNRAS, 303, L7 

Johnson, H. L. 1966, ARA&A, 4, 193 



-51 - 



Karachentsev, I. D., Karachentseva, V. E., Kudrya, Y. N., Sharina, M. E., & Parnovsky, S. 
L. 1999, BSAO, 47, 5 

Kautsch, S. J., Grebel, E. K., Barazza, F. D., & Gallagher, J. S. Ill, 2006, A&A, 445, 765 
Kennicutt, R. C. Jr. 1989, ApJ, 344, 685 
Kennicutt, R. C. Jr. 1998, ARA&A, 36, 189 

Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, MNRAS, 363, 2 

Kewley L. J., Geller, M. J., Jansen, R. A., & Dopita, M. A. 2002, ApJ, 124, 3135 

Kim, W.-T. & Ostriker, E. C. 2006, ApJ, 646, 213 

Kormendy, J. & Fisher, D. B. 2005, RevMexAA, 23, 101 

Landolt, A. U. 1992, AJ, 104, 340 

Levine, S. E. & Sparke, L. S. 1998, ApJ, 503, 125 

Lockman, F. J. 1977, AJ, 82, 408 

Mac Low, M.-M. & McCray R. 1988, ApJ, 324, 776 

Martos, M. A. & Cox, D. R 1998, ApJ, 509, 703 

Masset, F. S. & Bureau, M. 2003, ApJ, 586, 152 

Masset, F. & Tagger, M. 1997, A&A, 318, 747 

Matthews, L. D. 2000, AJ, 120, 1764 

Matthews, L. D. & Gallagher, J. S. 111. 1997, AJ, 114, 1899 
Matthews, L. D., Gallagher, J. S. Ill, & van Driel, W. 1999, AJ, 118, 2751 
Matthews, L. D., Gao, Y., Uson, J. M., & Combes, F. 2005, AJ, 129, 1849 
Matthews, L. D., van Driel, W., & Gallagher, J. S. Ill, 1998, AJ, 116, 1169 
Matthews, L. D. & Wood, K. 2003, ApJ, 593, 721 

Meier, D. D., Turner, J. L., Crosthwaite, L. P., & Beck, S. C. 2001, ApJ, 121, 740 
Meurer, G. R., Carignan, C, Beaulieu, S. F., & Freeman, K. C. 1996, AJ, 111, 1551 
Miller, S. T. & Veilleux, S. 2003, ApJ, 592, 79 
Nelson, A. H. 1976, MNRAS, 174, 661 
Noguchi, M. 2001, ApJ, 555, 289 

Noordermeer, E., Sparke, L. S., & Levine, S. E. 2001, MNRAS, 328, 1064 
Oiling, R. P. 1996, AJ, 112, 457 



-52- 



Oosterloo, T. 1993, A&A, 272, 389 

Papaderos, P., Loose, H.-H., Pricke, K. J., & Thuan, T. X. 1996, A&A, 314, 59 

Pietrznski, G., Ulaczyk, K., Gieren, W., Bresolin, P., & Kudritzki, R. P. 2005, A&A, 440, 
783 

Pohlen, M., Dettmar, R.-J., Liitticke, R., & Schwarzkopf, U. 2000, A&AS, 144, 405 

Quinn, P. J. 1987 in Nearly Normal Galaxies: From the Planck Time to the Present, ed. S. 
M. Faber, (New York: Springer), 138 

Quiroga, R. J. 1974, Ap&SS, 27, 323 

Quiroga, R. J. 1977, Ap&SS, 50, 281 

Rand, R. J. 1996, ApJ, 462, 712 

Reshetnikov, V. & Combes, P. 1997, A&A, 324, 80 

Revaz, Y. & Pfenniger, D. 2004, A&A, 425, 67 

Richter, O.-G. & Sancisi, R. 1994, A&A, 290, L9 

Saha, K. & Jog, C. J. 2006, A&A, 446, 897 

Salzer, J. J. & Norton, S. A. 1999, in The Low Surface Brightness Universe, lAU GoUoq. 
171, ed. J. I. Davies, C. Impey, and S. Phillips, (ASP: San Francisco), 253 

Sanchez-Saavedra, M. L., Battaner, E., Guijarro, A., Lopez-Gorredoira, M., & Gastro- 
Rodriguez, N. 2003, A&A, 399, 457 

Sancisi, R. & Allen, R. J. 1979, A&A, 74, 73 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Schneider, S. E. & Salpeter, E. E. 1992, ApJ, 385, 32 

Schoenmakers, R. H. M., Franx, M., & de Zeeuw, P. T. 1997, MNRAS, 292, 349 

Schwarzkopf, U. & Dettmar, R.-J. 2000, A&A, 361, 451 

Sellwood, J. A., Nelson, R. W., & Tremaine, S. 1998, ApJ, 506, 590 

Seth, A. G., Dalcanton, J. J., & De Jong, R. S. 2005, AJ, 129, 1331 

Seth, A. G., Dalcanton, J. J., Hodge, P. W., & Debattista, V. P. 2006, AJ, 132, 2539 

Sharina, M. E., Karachentsev, 1. D., & Tikhonov, N. A. 1999, Astr. Lett., 25, 322 

Shen, J. & Sellwood, J. A. 2006, MNRAS, 370, 2 

Sparke, L. S. 1995, ApJ, 439, 42 

Spicker, J. & Feitzinger, J. V. 1986, A&A, 163, 43 



- 53 - 



Stil, J. & Israel, F. P. 2002a, A&A, 389, 42 

Stil, J. & Israel, F. P. 2002b, A&A, 392, 473 

Swaters, R. A. & Balcells, M. 2002, A&A, 390, 863 

Swaters, R. A., Sancisi, R., & van der Hulst, J. M. 1997, ApJ, 491, 140 

Swaters, R. A., Schoenmakers, R. H. M., Sancisi, R., & van Albada, T. S. 1999, MNRAS, 
304, 330 

Swaters, R. A., van Albada, T. S., van der Hulst, J. M., & Sancisi, R. 2002, A&A, 390, 829 
Taylor, C. L. 1997, ApJ, 480, 524 

Taylor, C, Brinks, E., & Skillman, E. D. 1993, AJ, 105, 128 

Telles, E. & Terlevich, R. 1995, MNRAS, 275, 1 

Thuan, T. X., Hibbard, J. E., & Levrier, F. 2004, AJ, 128, 617 

Thuan, T. X. & Martin, G. E. 1981, ApJ, 247, 823 

Tifft, W. G. 1990, ApJS, 73, 603 

Tikhonov, N. A., Galazutdinova, O. A., & Drozdovsky, I. O. 2006, A&A, submitted ( |astro-pli/ 0603457D 

Toomre, A. 1981, in The Structure and Evolution of Normal Galaxies, Proc. Adv. Study 
Inst., (Cambridge University Press: Cambridge), 111 

Toomre, A. & Toomre, J. 1972, ApJ, 178, 623 

Tully, R. B. 1988, Nearby Galaxies Catalog (Cambridge University Press: Cambridge) 

TuUy, R. B., Pierce, M. J., Huang, J.-S., Saunders, W., Verheijen, M. A. W., & Witchalls, 
P. L. 1998, ApJ, 115, 2264 

Uson, J. M. 2007, VIA Test Memo #237 

Uson, J. M. & Matthews, L. D. 2003, AJ, 125, 2455 (Paper I) 

van der Hulst, J. M., Skillman, E. D., Smith, T. R., Bothun, G. D., McGaugh, S. S., & de 
Blok, W. J. G. 1993, AJ, 106, 548 

van der Kruit, P. C. & Shostak, G. S. 1984, A&A, 134, 258 

van Moorsel, G. A. & Wells, D. C. 1985, AJ, 90, 1038 

van Zee, L., Salzer, J. J., & Skillman, E. D. 2001, AJ, 122, 121 

van Zee, L., Skillman, E.D., & Salzer, J. J. 1998, AJ, 116, 1186 

Velazquez, H. & White, S. D. M. 1999, MNRAS, 304, 254 

Wada, K., Meurer, G., & Norman, C. A. 2002, ApJ, 577 



-54- 



Walker, I. R., Mihos, J. C, & Hernquist, L. 1996, ApJ, 460, 121 
Wang, J., Heckman, T. M., & Lehnert, M. D. 1999, ApJ, 515, 97 
Weinberg, M. D. 1995, ApJ, 455, L31 

Weinreb, S., Balister, M., Maas, S., & Napier, P. J. 1977, IEEE Trans. Microwave Theory 
Tech., vol. MTT-25, 243 

Wilcots, E. M. & Miller, B. W. 1998, AJ, 116, 2363 

Wilcots, E. M. & Prescott, M. K. M. 2004, AJ, 127, 1900 

Wu, Y., Charmandaris, V., Hao, L., Brandl, B. R., Bernard-Salas, J., Spoon, H. W. W., & 
Houck, J. R. 2006, ApJ, 639, 157 

Yuan, Q.-R. & Zhu, C.-X. 2004, ChA&A, 28, 127 

Zaritsky, D. & Rix, H.-W. 1997, ApJ, 477, 118 



This preprint was prepared with the A AS IM^jX macros v5.2. 



46°00' 



55' 



o 
o 
o 

CM 



g 



i 50' 





45' 



40" 



0.0499 



6:0503 



r ■ 



0.0458 



\ 



MGC2S37A 
0.0400 



; 0:2873 ; 



NGC 2537 
0.00 IB 



0.0479 
/ 



t— 



:0.D494 



V # 

0^605 




2973 



0.0540 
/ 



02484 



J:. 



0,0489 



45^ §£l^ 
Right Ascension (J2000) 



Fig. 1. — Image of the IC2233/NGC2537 field obtained from a mosaic of archival data from 
the SDSS. The grey-scale representation is based on photometric data from the u, g and r 
bands. The overplotted redshifts were also obtained from the SDSS archive. 
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Fig. 2. — WIYN images of IC 2233 in the i?-band (left) and continuum-subtracted Hq;+[Nii] 
(right). The small white rectangle on the left image denotes the kinematic center of the 
galaxy (see §6.1.1i) . Bleeding CCD columns from a bright foreground star are visible along 
the north-south direction near ajsooo = 08''13'"59^ (5j2ooo = +45°45'30". 
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Fig. 3. — Image of the IC 2233/NGC 2537 field in the 21-cm continuum, showing the central 
20' X 20' region. The position of IC 2233 is indicated by an ellipse and that of NGC 2537 by 
a circle. The contours are (-1.4 [absent],-!, 1,1.4,2,2. 8,4,. . . ,1024,1408)x0.46 mJy beam"^, 
where the lowest contour is 4o". We detect continuum emission from NGC 2537 but no 
significant emission from IC 2233. 
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Fig. 4. — Hi channel images for IC 2233. The spatial resolution is ~ 16". The contours 
shown are (— 2[absent], — 1, 1, 2, ...32) x 1.0 mJy beam~^. Only the channels showing line 
emission, plus one additional channel on either side, are plotted. Each panel is labeled with 
its corresponding heliocentric velocity. The systemic velocity is 553.4 ± 1.0 km s~^. The 
kinematic center of the galaxy ( §6.1.11) is indicated in each panel by a cross. 
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Fig. 5. — Results of two-dimensional elliptical Gaussian fits to the Hi emission of IC 2233 in 
several channels. Full-width at half-maximum of the major axis (after deconvolution of the 
synthesized beam) is plotted as a function of channel. The best fit parabola to the results 
for channels 30 - 37 is overplotted. The fit yields a minimum of the width for "channel" 
34.25, corresponding to a systemic velocity of 553.4 ± 1.0 km s~^ (see §6.21 for details). 
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Fig. 6. — Global (spatially integrated) H I profile of IC 2233 derived from our VLA data. 
The la error bars account for the statistical noise in each channel but do not include the 1% 
global calibration uncertainty. The inset shows the profile folded about the systemic velocity 
in order to highlight the difference in slope of the two edges. 
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Fig. 7. — H I total intensity contours of IC 2233 overlaid on the Ha+[N ll] image of the galaxy 
from Figure [2J Hi contour levels are (1, 1.4, 2, 2.8, 4. ..64, 78) x 18 mJy beam^^ m s~^. The 
highest contour departs from the "\/2" sequence in order to highlight the correspondence 
between local maxima in the H I column density and the location of Ha emission features. 
The synthesized beam is nearly circular with a FWHM ~ 16". 
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Fig. 8. — Deprojected Hi surface density profile for IC 2233, computed using an Abel 
inversion technique. H I surface density in units of solar masses per square parsec is plotted 
as a function of projected radius in kpc. The southern (approaching) side of the galaxy is 
shown as a dotted line and the northern (receding) side is shown as solid line. Values at very 
small radii (r <1.5 kpc) are highly uncertain (see §6.3p . 
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Fig. 9. — H I warp curve for IC 2233 derived by fitting Gaussians to a series of H I intensity 
profiles perpendicular to the major axis of the disk. A major axis positional angle of 172° 
was used. The displacement of the fitted centroid from the mean disk plane in arcseconds 
is plotted as a function of projected radius in arcseconds. The la error bars indicate the 
formal uncertainty in the fit. Arrows along the bottom of the figure denote the observed 
extent of the stellar disk. 
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Fig. 10. — Hi velocity field for IC 2233. Isovelocity contours ranging from 441-668 km 
are shown at 10.4 km increments (i.e., every other channel), overplotted on a linear 
grayscale representation of the same map. The greyscale range is 400-700 km s~^, with 
darker colors corresponding to higher velocities. 
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Fig. 11. — Hi position- velocity plots parallel to the major axis of IC 2233. The cuts were 
extracted in 3" strips along the disk major axis (z — 0), and sX z = ±15". The contours are 
(-1.4[absent], -1,1, 1.4, 2, 2.8, 4, 5.6, 8, 11,16,22,32,44) x 0.8 mJy beam"!. The greyscale 
range is 0-45 mJy beam~^. 




Fig. 12. — Sample minor axis position- velocity plots extracted at 30" intervals along the ma- 
jor axis of IC 2233. The offset along the major axis relative to the galaxy center is indicated 
in the upper right of each panel. The horizontal axes show radial velocity relative to the 
systemic velocity, and the vertical axes show distance from the midplane. Positive values of z 
correspond to the east side of the galaxy. Contour levels are (-2[absent],-l,l,2,4,8,16,32) xO.8 
mJy beam~^. 
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Fig. 13. — Hi rotation curve for IC 2233 (filled circles). The overplotted crosses are data 
points derived by Goad & Roberts 1981 using Ha emission-line spectroscopy. The axes are 
projected distance along the major axis in arcseconds and rotational velocity in km s~^. The 
dashed line indicates the systemic velocity of the galaxy. The location of the observed edges 
of the stellar disk are indicated with arrows. 
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Fig. 14. — Hi channel images for NGC 2537. Contours are overplotted on a greyscale 
representation of the same data. Contour levels are (— 2[absent],-l,l,2,4) xO.9 mJy beam~^. 
The greyscale range is — 5.0 mJy beam~^. Only the channels showing line emission, plus one 
additional channel on either side, are plotted. Each panel is labeled with its corresponding 
heliocentric velocity. 



-70- 




Fig. 15. — Global (spatially integrated) Hi profile for NGC 2537 derived from VLA data. 
The la error bars correspond to the statistical noise in each channel but do not include the 
~1% global calibration uncertainty. 




Fig. 16. — Hi moment maps of NGC 2537. For all data shown, the synthe- 
sized beam is roughly circular with a FWHM ~ 16". (a) Hi total intensity con- 
tours overplotted on a greyscale representation of the same data. Contour levels are 
(0.24,0.34,0.48,0.68,0.96,1.36,1.92) X 100 Jy beam"^ m s^^ The greyscale is linear from 0- 
200 Jy beam~^ m s~^. (b) Same as (a), but the Hi contours are overlaid on a blue image 
of the galaxy from the DSS. (c) Hi velocity field. Isovelocity contours ranging from 368- 
502 km are shown at 10.4 km increments (i.e., every other channel), overplotted on a 
linear greyscale representation of the same map. The greyscale range is 380-550 km s^^, with 
darker colors corresponding to higher velocities, (d) Hi velocity dispersion. The greyscale 
range is 2-20 km s~^. Maximum observed values are ~25 km s~^. 
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Fig. 17. — Azimuthally-averaged Hi surface density profile for NGC 2537. Mean Hi surface 
density measured in a series of concentric, elliptical annuli (in units of solar masses per square 
parsec) is plotted as a function of radius (in kpc). The error bars are statistical only and 
take into account the number of synthesized beams in each annulus as well as the number 
of channels that contribute to the H I signal at each position. Because of the nearly face-on 
orientation of the galaxy, no correction has been applied for inclination. The dashed line 
shows the critical density for star formation according to the criterion of Kcnnicutt (1989; 
see text). The dotted line shows the observed density scaled to account for the mass of 
helium and molecular gas. An arrow denotes the edge of the stellar disk. 
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Fig. 18. — Hi P-V plot along the major axis of NGC 2537. The data were averaged along 
a 15"- wide strip. Contours are (— 2 [absent],— 1,1,2,4,6,8) x 0.6 mJy beam~^. The grey scale 
range is 0—4.5 mJy beain"^. 
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Fig. 19. — Hi rotation curve for NGC 2537. No corrections for inclination or asymmetric 
drift have been applied to the observed circular velocities. The dashed line indicates the 
galaxy systemic velocity derived from a tilted ring fit (§8.2.4). 



